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ABSTRACT 

In  this  thesis,  the  theory  of  associated  polynomials  is 
applied  to  linear  circuits  and  systems.   Associated  polyno- 
mials are  defined  and  used  in  the  presentation  of  a  computer 
algorithm  suitable  for  accomplishing  a  partial  fraction  ex- 
pansion including  repeated  roots.   Variations  of  the  coeffi- 
cients of  a  polynomial  are  related  to  the  variation  of  its 
roots  (and  vice-versa)  by  associated  polynomials.   These  re- 
sults are  used  in  a  root  solving  process  and  to  exoress  sen- 
sitivity coefficients  in;  a  new  analytical  form.   Lastly ,  a 
theory  of  compensating  parameter  adjustments  is  developed 
with  applications  to  self-adaptive  systems.   It  is  shown  that 
by  compensating  adjustments  of  any  two  parameters  in  the  sys- 
tem which  are  linearly  related  to  the  coefficients  of  the 
characteristic  equation,  adaptive  compensation  is  possible 
for  a  real  root  of  the  closed  loop  transfer  function  only 
and  that  for  complex  roots,  the  real  oT     imaginary  parts  may 
be  kept  invarient  but  not  both. 
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CHAPTER  I 
Associated  Polynomial  Theory 
1.1:  Introduction 

This  thesis  is  concerned  with  the  investigation  and 
application  of  associated  polynomials  to  the  analysis  and 
synthesis  of  linear  systems.   As  is  shown,  every  polyno- 
mial of  degree  n  can- be  used  to  generate  a  sequence  of  n+1 
polynomials  of  decreasing  degree  called  the  associated 
polynomials  of  the  defining  polynomial.   The  associated 
polynomials  depend  in  a  particularly  simple  wav  uoon  the 
coefficients  of  the  original  polynomial.   Due  to  the  rela- 
tionship between  the  coefficients  of  the  original  polyno- 
mial, the  associated  polynomials,  and  the  properties  of  the 
original  polynomial  which  can  be  expressed  simply  in  terms 
of  the  associated  polynomials,  several  useful  applications 
and  results  are  possible  in  network  analysis  and  control 
theory.  _, 

The  concept  and  use  of  associated  polynomials  is  not  new 
in  that  these  polynomials  have  appeared  in  many  guises  in  the 
literature    ,  usually  related  to  some  particular  application. 
Associated  polynomials  first  appeared  in  the  literature  in 

the  late  eighteenth  century  when  Legrange  used  them  to  obtain 

(2) 
a  solution  of  the  homogeneous  linear  difference  equation 

(1  3) 

Associated  polynomials  have  been  called  Horner  polynomials 

as  they  were  used  directly  in  Horner's  method  or  rule    .  pk& 
the  purely  mathematical  aspects  of  associated  polynomials  &te 
not  of  direct  significance  to  this  thesis.   The  use  of  asso- 
ciated polynomials  in  Horner's  method  is  presented  separately 
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in  appendix  A.   What  often  occurs  in  engineering  applications 
is  not  an  associated  polynomial  itself,  but  rather  a  number 
corresponding  to  an  associated  polynomial  evaluated  at  a 
particular  root  of  the  original  polynomial. 

Recently    ,  the  properties  of  associated  polynomials 
have  been  investigated  in  themselves  and  new  concepts  and  re- 
sults introduced.   Reference  1  is  written  in  a  pure  mathemat- 
ical language  with  only  mathematical  formulas  as  results. 
The  object  of  this  thesis  is  to  demonstrate  that  associated 
polynomials  have  very  useful  properties  and  applications  in 
the  fields  of  circuit  and  systems  theory. 

This  thesis,  divided  into  five  chapters,  presents  the 
basic  associated  polynomial  theory  in  chapter  one.   In  chapter 
two,  this  basic  theory  is  used  to  present  an  original  algo- 
rithm suitable  for  computer  programming  for  accomplishing  a 
partial  fraction  expansion  including  repeated  roots  of  high 
order  of  multiplicity.   Two  theorems  are  developed  in  chapter 
three  that  show  that  the  variation  of  the  coefficients  of  a 
polynomial  with  the  variation  of  its  roots  or  vice-versa  is 
determined  by  associated  polynomials.   The  theorems  are  made 
possible  by  deriving  the  orthnormality  relationship  between 
the  set  of  associated  polynomials  and  the  roots  of  the  defin- 
ing polynomial.   The  above  theorems  have  important  engineering 
applications  in  that  they  are  used  in  a  root  solving  process 
in  chapter  three,  and  in  chapter  four  to  express  sensitivity 
coefficients  in  a  new  analytical  form.   Chapter  five  shows 
how  associated  polynomial  theory  may  be  used  in  the  develop- 
ment of  a  theory  of  compensating  parameter  adjustments  with 
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application  to  self-adaptive  systems  in  which  a  change  in  a 
variable  parameter  in  a  plant  may  be  compensated  for  bv  a 
change  in  another  variable  parameter  such  that  the  dominant 
root  performance  of  the  overall  system  is  controlled.   An 
important  theorem  of  this  chapter  is  the  proof  that  if  a  feed- 
back control  system  has  a  real  root  then  the  location  of  this 
root  may  be  kept  constant  by  compensating  adjustments  of  any 
two  parameters  in  the  system  which  are  linearly  related  to 
the  coefficients  of  the  characteristic  equation.   It  is  shown 
that  when  the  compensating  parameters  are  associated  with 
simple  zeros  and  poles  of  the  open  loop  transfer  function, 
such  adaptive  compensation  is  possible  for  a  real  root  of  the 
closed  loop  transfer  function  only  and  that  for  complex  con- 
jugate roots  the  real  or  imaginary  parts  may  be  kept  invar- 
ient  but  not  both.   This  theorem  is  of  fundamental  importance 
in  adaptive  control  theory. 
1.2:  Definition  of  Associated  Polynomials 

Let  a  monic  polynomial,  F(s) ,  of  degree  n  be  expressed 
as  : 


-l  n 

F(s)  =  a  sn  +  a   , sn  X  +  . . .  +  an  =   J  a.s^  ;  a   =  1     (1-1) 
n      n-1  0    *  j  n 


The  condition  that  F(s)  be  monic,  i.e.,  a   =  1 ,  is  convenient 
rather  than  essential.   Let  the  distinct  roots  of  F(s)  be 
labeled,  p.,  j  =  l,2,...,n. 

If  u  is  any  variable  different  from  s,  then  the  polyno- 
mials A.(u);  j  =  0,1,..., n  are  defined  by  the  operation: 

sF(s)-uF(u)  _  _  ,.     .    .  ,  ,   .      ^_    -  ,  s  n  _   r  a  f„\«J 


=  An(u)  +  A, (u)s  +  ...  +  A  (u)sn  =   I  A.(u)s 
0        1  n         _> L  n    1 

(1-2) 
11 


s-u         0        1  n         >Q 


The  set  of  polynomials,  A,(u);  j  =  0,1,..., n  are  called  the 
associated  polynomials  of  F(s). 

Example  1-1 

2 
Given  the  polynomial:   F(s)  =  s   +3s+2 

The  associated  polynomials  of  F(s)  are  found  by  the  following 

operations : 

F(u)  =  u2  +  3u  +  2 

Performing  the  multiplication  and  division  of  equation  1-2, 

sF(s)-uF(u)     2    ,_   .     ,  2l0  l~ls 

! — - : — -  =  s   +  (3+u)s  +  (u  +3u+2) 

s-u 

The  associated  polynomials  of  F(s)  are: 

AQ(u)  =  u2  +  3u  +  2 

A1 (u)  =  3  +  u 

A2(u)  =  1 

The  basic  definition  of  associated  polynomials  given  bv 
equation  1-2  may  be  simplified  by  noting  the  relationship  of 
associated  polynomials  and  the  coefficients  of  the  basic 
polynomial.   This  relationship  is  given  by  Theorem  1-1. 
THEOREM  1-1 

—  +- Vi 

The  associated  polynomials  of  an  n     order  polynomial, 
F(s) ,  may  be  found  from  the  coefficients  of  F(s)  by  the  recur- 
sion formula: 


n 

A.(u)  =  I    a.u  ^     :  j  =  0,1,. ..,n 
3  i=j 
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Proof: 

Multiplying  equation  1-2  by  (s-u) ; 

sF(s)-uF(u)  =  [AA(u)+A,  (u)s+...+A  (u)sn](s-u)     (1-3) 

u     l  n 

Substituting  equation  1-1  for  F(s)  and  multiplying  throuqh 
by  (s-u),  equation  1-3  becomes: 

n+1        n  ■/.  i    «  /  \  n+1 

as     +  a   ,s   +  ...  +  ans  -  uF(u)  =  A  (u)s    + 
n        n-1  0  n 

[A^  , (u) -UA^ (u) ]s  +  ...  +  [An (u) -uA, (u) ]s  -  uAn(u) 
n-i      n  u      1  u 

(1-4) 

Comparing  the  coefficients  of  the  equal  powers  of  s, 

a„    =  A  (u)  (l-5a) 

n      n 

a„  ,  =  A^  , (u)  -  uA  (u)  (l-5b) 

n-l    n-l        n 

aB  -  =  A^  0(u)  -  uA   . (u)  (l-5c) 

n-z  n-z  n-l 


aQ   =  AQ(u)  -  uAx(u)  (l-5d) 

Since  a  =1/  equation  1-5  may  be  explicitly  expressed  as: 


A.(u)  =  uA.+1(u)  +  a.  :  j  =  0,l,...,n-l 


(1-6) 


By  definition  of  the  set  of  associated  polynomials,  A.(u),  j 
ranges  from  0  to  n.   Hence  let: 

A  .,  (u)  =  0  (1-7) 

n+1 

Using  equation  1-7,  equation  1-6  becomes: 

A.(u)  =  uA.+1(u)  +  a.  :  j  =  0,1,. . . ,n  (1-8) 
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In  order  to  form  an  explicit  expression  for  the  associated 

polynomials,  note: 

A  (u)    =  a^  (l-9a) 

n         n 

A^  -  (u)  =  uA  (u)  +  a^  ,  =  a   ,  +  a  u  (l-9b) 

n-l        n       n-1    n-1     n 

2 
A^  0(u)  =  uA   . (u)  +  a^  0  =  a^  0  +  aM  , u  +  a  u    (l-9c) 
n-2,  n-1       n-2     n-2    n-1     n 


Art  (u)    =  uA,  (u)  +  an  =  an   +  a,  u  +  ...  +  a  u      (l-9n) 
0  1        U     0     1  n 


for  j  =  0 


n 

l 


A  (u)  =  I    a.u  (l-10a) 

°       i=0  x 


for  j  =  1: 

n     .  , 
A,  (u)  =   J"  a.u1  X  (l-10b) 

hence : 

n 
A.(u)  =  I    a.u1  :    j  =  0,1,. ..,ri         (1-11) 

Corollary  1-1 

The  associated  polynomial,  A.(u) ,  is  a  polynomial  of 
degree  n  -  j . 
Corollary  1-2 

The  particular  associated  polynomial,  AA(u),  is  always: 

AQ(u)  =  F(u)  (1-12) 

Example  1-2 

Using  the  results  of  Theorem  1-1,  the  associated  polyno- 
mials of  the  polynomial,  F(s),  of  example  1-1  may  be  deter- 
mined as  follows: 
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F(s)  =  s   +  3s  +  2  =  a2s   +  a,s  +  aQ 

n     .  . 
A.  (u)  =  I    aj_u1  3     :  j  =  0,1,2 
J      i=j 

The  associated  polynomials  of  F(s)  are: 
A2(u)  =  a2  =  1 

A, (u)  =  a,  +  a2u  =  3  +  u 

2  2 

AQ(u)  =  aQ  +  a,u  +  a2u  =  2  +  3u  +  u 

which  agrees  with  example  1-1. 

1.3:  Characteristics  of  Associated  Polynomials 

The  above  section  has  shown  the  elementary  relationship 
of  associated  polynomials  to  the  coefficients  of  the  original 
polynomial.   There  are  other  interesting  properties  of  asso- 
ciated polynomials  that  can  now  be  considered.   These  prop- 
erties include  the  relationship  between  associated  polyno- 
mials and  the  original  polynomial  which  leads  to  the  concept 
of  synthetic  division  and  the  fact  that  associated  polyno- 
mials are  elementary  symmetric  functions.   These  properties 
are  stated  in  the  following  theorems. 
Theorem  1-2 

The  relationship  between  a  polynomial/  F(s),  and  its 
associated  polynomials,  given  any  particular  root,  p.,  of 
F  ( s )  i  s  : 

Zi*L=    ?a-(p-)sj_1 
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Proof 


For  convenience,  equation  1-2  is  repeated  as: 


n 


(sF(s)) [s,u]    =   A0(u)    +   Ax(u)s   +    ...    +   An(u)s  (1-2) 


n-1. 


Equation    1-2   may   be    rewritten   as: 
(sF(s))  [s,u]    =   AQ(u)    +    s[A,(u)    +   A2(u)s   +    ...    +   A    (ujs11-"1"] 

(1-13) 
However  the  expression  in  the  brackets  of  equation  1-13  may 
be    equated   to: 

F(S)"F(U)    =   A,(u)    +   A9(u)s    +    ...    +   A    (u)s11"1    =      I    A.(u)sj_1 
s-u  1  z  n  ■    i .  i 

1  =  1    J 

(1-14) 
Since  F(u)  =  A~(u) ,  equation  1-14  may  be  rewritten  as: 


F(s) 
s-u 


AQ(u)     n 


■ +  I   A.(u) 

i  =  l  ^ 


,j-l 


(1-15) 


Note  that  equation  1-15  states  that  the  quotient  and  remainder 
after  division  by  a  monic  linear  polynomial  may  be  calculated 
recursively.   As  shown  in  appendix  A,  the  nature  of  this  re- 
cursion leads  to  synthetic  division. 

Let  u  =  p. ,  be  a  root  of  the  polynomial,  F(s) ,  then, 


since  F (p . )  =  0 : 


AQ(u) 


s-u 


F(u) 


u=p. 


u=p.  s-u 

and  equation  1-15  may  be  rewritten  as : 

n 


F(Pj>  =  0 
s-p. 


(1-16) 


F(s)  _ 
s 


ii^-   I  A.  (p. 

"Pi   j=i  5  pi 


)s 


j-1 


(1-17) 


Theorem  1-2  may  be  numerically  verified  by  the  follow- 
ing example. 
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Example  1-3 
Given: 

F(s)  =  s4  +  10s3  +  35s2  +  50s  +  24  =  (s+1)  (s+2)  (s+3)  (s+4) 
Let: 

P1   =  _1?  P2  =  ~2;  P3  =  ~3'*  P4  =  ~4* 
The  coefficients  of  F(s)  are: 

aQ  =  24;  a,  =  50;  a2  =  35;  a3  =  10;  a4  =  1. 

The  associated  polynomials  of  F(s)  are: 
A4(u)  =  an  =  1 

A3(u)  =  a3  +  a.u  =  10  +  u 

2  2 

A2(u)  =  a2  +  a3u  +  a.u  =  35  +  lOu  +  u 

2       3  2     3 

A, (u)  =  a,  +  a2u  +  a3u  +  a.u  =  50  +  35u  +  lOu  +  u 

Note:   The  particular  trivial  associated  polynomial  AQ(u)  is 
no  longer  required  and  is  not  calculated  in  future  examples. 
Consider  the  root,  p2 ,  of  F(s),  then: 

A1(p2)  m    12;  A2(p2)  =  19;  A3(p2)  =  8;  A4(p2)  =  1. 

Expanding  equation  1-17  and  evaluating: 

|||i  =  Al(p2)  +   A2(p2)s  +  A3(p2)s2  +  A4(p2)s3 

=  12  +  19s  +  8s2  +  s3 


s+2 


As  a  check : 


(s3+8s2+19s+12)  (s+2)  =  s   +  10s3  +  35s2  +  50s  +  24   as  given 
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Corollary  1-3 

Theorem  1-2  developed  a  relationship  between  a  polynomial, 
F(s),  and  its  associated  polynomials  given  a  Darticular  root 
of  F(s).   This  relationship  may  be  extended  to  polynomials 
with  multiple  roots  so  that  the  following  formulas  are  avail- 
able. 
CASE  (1) :   TWO  REPEATED  ROOTS 

'  -,  =  A  (p.)sn"2  +  [A       ,  (p.)+p.A  (p.Hs11"3  + 
2     n  c  l  n-l   l   -  i  n  -  l 


(s-p.) 


[A   0(p.)+D.A   ,(p.)+p.2A  (p.)]sn~4  +  ... 
n-2  *i    i  n-l  *i   ri   n  --i 

+  [A9(p.)+p.A.,(p.)  +  ...+p.n~2A  (p.)]      (1-18) 


Expressed  as  an  explicit  relationship: 

F(s)      nv2   h    ?     j-(h+2)_  ,   ,       ,•-  1Q> 

2=  I      s     2-   Pi        A.  (p.)       (1-19) 

(s-p.)     h=0     j=h+2  1        :   1 

CASE     (2) :       THREE    REPEATED    ROOTS 

F(s)  -     /       v    n-3         ri  /       \  ,  o      TV     /       \  i    n-4 

-=■  =   A    (p.  )  s  +    [A      ,  (p. )+2p. A    (p.)]s 

3  n   *i  n-l   ^i        *i    n    -i 


(s-Pi) 

+     [A      0(p.)+2p.A      . (p. )+3(p. )2A    (p.)]sn"5+    ... 
n-2  f l         -i   n-l   ^i  *i        n   ^i 

+     [A3(pi)+2piA4(pi)+.. .+(n-2) (pi)n"3An(pi)] 

(1-20) 
Expressed  as  an  explicit  relationship 

n-2   ,  ,    n 


F(s) 


=   I   s*"1    I       [j-(h+l)]p.^(h+2)A  (p.)   (1-21) 


(s-p.)3    h=l      j=h+2 
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CASE  (3) :   FOUR  REPEATED  ROOTS 

-— T-  VPi>sn-4   +    [An.1(Pi)  +  3oiAn(Di)]sn-5 
(s-Pi) 

+    tAn-2(Pi)+3PiVl(Pi)+6(Pi)2An(Pi)lsn"6 

+[An-3(Pi)+3PiAn-2(Pi)+6(Pi)2An-l(Pi)+10(Di)3An(Pi)lsn"7 

+    ...    +    [A.(p,)+3p,At.(p.)  +  ...  +  [l+...  +  (n-3)]p.n~4A    (o.)  ] 
41  x  5     1  ~i  ni 

(1-22) 

Expressed  as  an  explicit  relationship, 

F(S)  A    =  T  a*"1  I       {l+...+  [i-(h+2)]}p.^(h+3)A.(p.) 
(s-p.)     h=l     j=h+3  L        J      1 

(1-23) 
Example  1-4 
Given: 

F(s)  =  s5  +  13s4  +  67s3  +  171s2  +  216s  +  108  =  (s+2)2(s+3)3 
Let:   p,  =  -2;  p2  =  -3. 

The  associated  polynomials  of  F(s)  are: 
A5(u)  =  1 

A4(u)  =  13  +  u 

A3(u)  =  67  +  13u  +  u2 

A2(u)  =  171  +  67u  +  13u2  +  u3 

A1(u)  =  216  +  171u  +  67u2  +  13u3  +  u4 

Evaluating  the  associated  polynomials  at  p, : 

A5(p1)  =  1;  A4(Pl)  =  11;  A3(P;L)  =  45;  A^p^  =  81. 
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Evaluating  the  associated  polynomials  at  p2  : 
A5(p2)  =  1;  A4(P2)  =  10;  A3(d2)  =  37. 

Applying  equation  1-18: 

F(S)  2  =  A5(p1)s3  +  [A4(p1)+p1A5(p1) ]s2  +  [A3(p1)+p1A4(p1)+ 
(s-px) 

2  2  3 

pl  A5(Pi^s  +  [A2(Pi)+PiA3(Pi)+Pi  A4(P1)+P1  A5(PX) 


Evaluating : 


Is3   +    [ll-2]s2   +    [45-22+4]s   +    [81-90+44-8] 


i  \2 


or: 


=  s3  +  9s2  +  27s  +  27  =  (s+3)3 


(s-p1) 


Applying  equation  1-21: 

(S)  3  =  A5(p2)s2  +  [A4(p2)+2p2A5(p2) ]s  +  [A3(p2)+2p2A4(p2) 


(s-p2) 


2 

+3p2  A5(p2) ] 


Evaluating : 


=    Is2    +    [10-6]s    +    [37-60+27] 


/     \  3 
(s-p2) 


F(S)    =  s2  +  4s  +  4  =  (s+2)2 


(s-p2) 


Theorem  1-3 


Any  result  which  involves  an  associated  polynomial, 
A.(p^),  evaluated  at  a  particular  root,  can  always  be  inter- 
preted as  a  result  involving  elementary  symmetric  functions, 
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i.e. 


Aj(pi)  =  (-l)jan_j  :  j  =  l,2,...,n-l 


—  4-Vi 

where  a  _ .  denotes  the  n- j  "   elementary  symmetric  function  of 
the  roots  of  F(s)  with  p.  absent. 

Proof: 

(5) 
Elementary  symmetric  functions  are  defined     as  fol- 
lows:  If  f(x)  is  a  polynomial  of  degree  n  with  complex  coef- 
ficients, the  functions  obtained  by  taking  the  sum  of  all  the 
roots  of  f (x) ,  the  sum  of  all  products  of  pairs  of  the  roots 
of  f (x) ,  the  sum  of  all  products  of  triplets  of  the  roots  of 
f(x)  and  so  on  to  the  product  of  all  the  roots  of  f(x)  may 
be  expressed  rationally  in  terms  of  the  coefficients  of  f(x). 
These  functions  of  the  roots  of  a  polynomial,  f(x),  are 
galled  elementary  symmetric  functions. 

Hence,  for  a  polynomial,  F(s),  of  degree  n  with  roots, 

]?1  /  P2  '  •  •  •  *  Pn  • 

,->  /  \    /     \  /     \    /     \     n      n— 1      n-2     ,  /  ,  s  n 
F(s)  =  (s-p, ) (s-p2) • • • (s-p  )  =  s   -  a,s    +  o 2s        -...+(-1)  a 


where 


(1-24) 

a,  =  p,  +  p~  +  ...  +  p  (l-25a) 

1    ^1    c 2  n 

°2   =   PlP2  +  P1P3  +  •*;+  pn-lpn         (l-25b) 


an  =  P1P2P3*  * ,pn  (l-25c) 

a, ,  0^,  ...,  a   are  the  elementary  symmetric  function  of  the 
roots,  p^,  p2,  . . . ,pn« 
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Rewriting  equation  1-17  as: 

G(s)  =  P&   =A1(p.)  +A,(p,)s+  ...  +  A  (p.)sn_1  (1-26) 
s-p ,     l^i      2.    c  l  n  -  l 

G(s)  is  a  polynomial  of  degree  n-1.  By  comparison  with  equa- 
tion 1-24,  the  associated  polynomials  of  F(s),  evaluated  at 
a  particular  root,  p. ,  can  be  seen  to  be  elementary  symmetric 
functions  of  p. ,  p~,  . ..,  p  with  the  particular  root,  p.  =  0 
Changing  the  notation  to  conform  with  the  already  estab- 
lished associated  polynomial  notation: 

An-l(pi*  =  "(Pl+P2+* * *+pn) ;  pi  =  °  (l-27a) 

An_2(pi)  =  + (p1p2+p1p3+. . -+Pn_iPn) 1    Pi  =  °     (l-27b) 

•    •    # 

Expressing  equations  1-27  explicitly: 


A. (p.)  =  (-l)^a    with:   p.  =  0  and  j  =  l,2,...,n-l 
j   c i  n-j        ri         J 


(1-28) 
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CHAPTER  II 
Partial  Fraction  Expansion  Using  Associated  Polynomials 
2.1:  Introduction 

The  basic  associated  polynomial  theory  developed  in  the 
preceding  chapter  has  immediate  application  in  the  field  of 
network  analysis  and  network  synthesis.   Specifically,  asso- 
ciated polynomials  may  be  used  as  an  alternate  method  for 
performing  the  Heaviside  or  partial  fraction  expansion  of  a 
ratio  of  polynomials. 

In  network  analysis,  a  method  of  finding  the  impulse  re- 
sponse, h(t),  from  the  system's  function  is  to  perform  a  par- 
tial fraction  expansion  of  the  response's  Laplace  transform, 
H(s),  and  then  use  standard  inverse  transform  tables.   Many 
network  synthesis  procedures  require  that  network  functions 
be  broken  up  into  the  sum  of  a  number  of  terms  which  can  be 
identified  with  network  elements.   In  the  case  of  the  well 
known  Cauer  procedures,  this  involves  partial  fraction  ex- 
pansions . 

The  use  of  associated  polynomials  in  partial  fraction 
expansions  is  not  the  easiest  method  available  when  hand  cal- 
culations are  to  be  used.   However,  the  calculations  involved 
using  associated  polynomials  are  very  adaptable  to  the  dig- 
ital computer.   Although  an  actual  computer  program  is  not 
presented  here,  this  chapter  presents  the  algorithm,  or  the 
step-by-step  procedure  for  accomplishing  a  partial  fraction 
expansion  of  any  ratio  of  polynomials. 
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Consider  the  rational  function: 


m 

1 
s 


c(s)  "  "n 7  "  5TiT  (2_1) 

I  ais± 

i=0  x 
where  a   =  1  and  n  _>  m.   If  equation  2-1  represents  the  trans- 

4- In 

formed  solution  of  an  n    order  differential  equation,  then 

the  denominator  term  is  the  characteristic  equation  of  the 

( 6 ) 
system  under  study.   It  can  be  shown     that  the  roots  of  the 

characteristic  equation  can  occur  in  no  more  complicated  form 

than  a  pair  of  complex  conjugate  roots  or  as  multiple  roots. 

Therefore,  for  convenience,  partial  fraction  expansions  using 

associated  polynomials  are  developed  using  three  different 

cases.   These  are: 

Case  (1):   C(s)  contain  first-order  poles  only, 

Case  (2):   C(s)  contains  a  pair  of  complex  poles  plus 
first-order  poles, 

Case  (3):   C(s)  contains  multiple  poles  plus  first  order 
poles . 
2.2:  Case  (1)  Partial  Fraction  Expansion 

If  the  roots  of  the  denominator  of  equation  2-1  are  each 
real  and  distinct,  then  by  the  fundamental  theorem  of  algebra, 
the  function,  C(s),  may  be  rewritten  as: 

C(S)  =  t T7 KTTT 7 r  (2-2) 

(s-p1) (s-p2) . . . (s-Pn) 
The  partial  fraction  expansion  of  equation  2-2  is: 

C(s)  =  -i ^-r-  +  t ir-  +  ...  +  -, ^t-        (2-3) 

(s-pjT    (s-p^  (s-Pn) 
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The  coefficient,  k.,  called  the  residue  of  C(s)  at  s  =  p . , 

i  =  1,2,... ^n,  is  an  undetermined  coefficient  of  the  expansion 

and  is  found  by: 

N(s) 


k.  = 

1 


DTiT 


(s-p. ) 


(2-4) 


s=p. 


Equation  2-4  is  the  classical  method  of  finding  the 
residues  of  a  partial  fraction  expansion  when  all  the  roots 
are  simple  and  distinct.   This  calculation  is  easy  to  per- 
form by  hand  calculations  but  requires  several  extra  steps 
on  a  digital  computer  which  will  not  divide  by  zero.   Asso- 
ciated polynomial  theory  may  be  used  instead  of  the  above 
computations  by  the  following  method. 

Rearrange  equation  2-4  as : 


k.  = 

l 


N(s) 
D(s)/(s-pi) 


i  =  1 ,2  , .  .  .  ,n 


(2-5) 


S-P. 


The  denominator  of  equation  2-5  can  be  expressed  using  equa- 
tion 1-17  as : 


D(s) 


n 


3-1; 


s-p 


=  I       S^XA,(p,) 


i   j  =  l 


1    "  i 


(2-6) 


where  A.  (p.)  is  the  set  of  associated  polynomials  of  D(s) 
evaluated  at  s  =  p..   Hence  equation  2-5  can  be  expressed  as: 


k.  = 

i 


N(Pi) 


n 


j-l; 


i  =  1 ,2 , . . . ,n     (2-7) 


I      PiD'XA  (p. 

i=l  1    :  1 


Example  2-1 

2 

1-.     j    ni    \         12s   +  22s  +  6   .        .  .  ,  j-    .  . 

Expand  C(s)  =  ^ ^ "Y  a  partial  fraction  expansion. 

s   +  3s   +  2s 
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Solution : 

rrc\    N(s)    12s2+22s+6     ^1        ^2__        k3  ' 

{S)         DTiT  "     3_,_,  2^0     s     s+1    s+2 
s  +  3s  +zs 

Let:   p,  =  0;  p„    =  -1;  p   =  -2. 

Form  the  associated  polynomials  of  the  denominator,  D(s): 

A3(u)  =  1 

A2(u)  =  3  +  u 

Ax(u)  =  2  +  3u  +  u2 

Applying  equation  2-7,  the  partial  fraction  coefficients  are 

.      N(s=0)     6    - 
*1  "  A1(0)+0     2    J 

k   =  N(s=  -1) 


where 


and: 


Hence : 


2    A1(-l)+(-l)A2(-l)+(-l)2A3(-l) 


A1(-l)  =  0;  A2(-l)  =  2;  A3(-l)  =  1 


N(s=  -1)  =  12  -  22  +  6  =  -4 


k   =   ~4    =  4 
K2    0-2+1    * 


k   =  N(s=  -2) 


where : 


and 


3    A1(-2)+(-2)A2(-2)+(-2)2A3(-2) 


A1(-2)  =  0;  A2(-2)  =  1;  A3(-2)  =  1 


N(s  =  -2)  =  48  -  44  +  6  =  10 
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Hence: 

.     10 _  - 

K3  ~  0+l(-2)+l(4)  ~  b 

The  partial  fraction  expansion  of  the  function,  C(s),  is 
nl    .    12s2+22s+6    3     4      5 

C(S)   =   — o = =   —  +   — -r-     + 


3^-,  2^„      s    s+1    s+2 
s  +3s  +2s 

2.3:  Case  (2)  Partial  Fraction  Expansion 

When  the  denominator  of  equation  2-1  contains  a  pair  of 

complex  roots  plus  first  order  poles,  then  the  function,  C(s), 

is  written: 

C(S)  =  (s-Pl) (s-Pl*) (s-p2) . . . (s-pn_x)        (2"8) 
The  partial  fraction  expansion  of  equation  2-8  is: 

JC^               ^*T                 Jv#*                            JC     ^ 
C(S)   =   I =-T  +   -7 ^-xy-  +  -, =-y-  +   ...   +  -> ^4r-   (2-9) 

(s-px)   Ts-pj*T    (s-p2)         (s-pn_1) 

As  shown  in  reference  6,  the  coefficients,  k.  and  k,*, 
must  be  complex  conjugates  and  either  may  be  found  by  the 
same  technique  used  in  case  (1)  above. 

Example  2-2 

4 
Expand  C(s)  =  * by  a  partial  fraction  expansion. 

s(s  +2s+2) 

Solution: 

c(a,  _  N(s)  _       4      _ .  £l    K_2 K2 

-  DTiT  "  s(s2+2s+2)  -'  s    s+l-j   STI7T 

Let  px  =  0;  p2  =  -1+j;  p2*  =  -1-j. 

Form  the  associated  polynomials  of  the  denominator,  D(s): 

A3(u)  =  1 
A2  (u)  =  2  +  u 
A3(u)  =  2  +  2u  +  u2 

27 


Applying  equation  2-7,  the  partial  fraction  coefficients  are 

_     N(0)  _  4 
Kl  "  A1(0)+0  ~  2  "  * 


k„  = 


N(-l+j) 


A1(-l+j)+(-l+j)A2(-l+j)+(-l+j)2A3(-l+j) 

Evaluating: 

A1(-l+j)  =2-2+2j-2j=0 

A2(-l+j)  =  1  +  j 

A3(-l+j)  =  1 

N(-l+j)   =  4 
hence : 

k2  "  ITT  =  -1  +  j 
As  the  coefficient  k2*  must  be  the  conjugate  of  k2  , 

v  ■  -±i 

The  partial  fraction  expansion  of  the  function,  C(s),  is: 

C(s)  =  i =  1  +   ~1+j  +   "1-J 

s(s2+20+2)    s    s+1^    s+1+^ 

2.4:  Case  (3)  Partial  Fraction  Expansion 

When  the  function,  C(s),  is  of  the  form: 

P/0x  _  N(s) _  N(s)       (0    ,m 

c(s)  -  - — rr; — ; — ; r _  ^r       {     ] 

(s-p1)  (s-p2) . . . (s-pn_R) 
The  partial  fraction  expansion  is: 

„,  .       m         m-1  1       2  n-K 

C  (s)  =  7T  +    77— s-  +  .  .  .  +  +  +  .  .  .  + 


(s-Pl)K  S"nn-K 

(2-11) 


(s-p1)     (s-p^^) 
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As  shown  by  equation  2-11  the  partial  fraction  expansion  is 
modified  to  include  as  many  additional  coefficients  as  the 

order  of  the  multiple  root. 

(6  1) 
Several  methods   '    are  available  for  the  evaluation 

of  the  coefficients  of  the  multiple  roots.   Computer  programs 

(8  9 ) 
are  also  available   '    that  perform  partial  fraction  expan- 
sions of  expressions  containing  multiple  roots  but  as  shown 
in  reference  7 ,  there  are  several  inherent  sources  of  error 
that  tend  to  degrade  the  numerical  accuracy  as  the  order  of 
multiplicity  of  the  root  increases. 

This  section  introduces  a  procedure  for  expanding  a  mul- 
tiple root  expansion,  such  as  equation  2-10,  using  a  combina- 
tion of  associated  polynomial  theory  and  simple  algebraic 
manipulations.   The  only  source  of  error  in  the  method  is 
that  the  denominator  polynomial,  D(s),  must  be  in  factored 
form. 

The  classical  method     of  determining  the  partial  frac- 
tion expansion  of  the  function,  C(s)  is: 

N    N. 

C(s)  =  I         f         °iK  v  (2-12) 

i=l  K=l  (s-p. ) 

where  the  coefficients,  cK,  are  given  by: 


-K         n. 

1    L  ->  s=p. 


1    d  i 


\LK    (N.-K)  !  ,  N.- 
l      ds  l 

■  l 

The  interdependence  of  the  successive  derivatives  in 

equation  2-13  is  a  possible  cause  of  errors  which  cumulate 

as  the  multiplicity,  N.,  increases.   A  second  method  of 
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determining  the  coefficients,  c . v ,  using  algebraic  manipula- 
tions  only  is  demonstrated  in  examples  2-3  and  2-4. 
Example  2-3 
Given  the  rational  function: 


C(s)  = 


Ti+ry 


(2-14) 


s(s+l)  (s+1) 

Using  the  classical  method  of  determining  the  coefficients, 
k,  and  k2,  are  found  to  be: 


k,  =  s  C(s) 


s=0 


k2  =  (s+ir  C(s) 


=  1 


s=-l 


H  1/-1  =  -1 


To  find  the  coefficient,  k3,  multiply  both  sides  of  equation 
2-14  by  the  denominator  of  k3,  i.e.,  (s+1).   Using  the  known 
values  of  k.  and  k2 ,  combine  the  terms  that  do  not  go  to  zero 
as  s-*— 1.   There  will  be  a  cancellation  of  the  undeterminate 
terms  after  the  combination  of  terms.   After  the  cancellation, 


let  s  =  -1  and  solve  for  k^. 


(s+1) 


s(s+l) 


=  (s+l)-f  + 


(s+l)k2    (s+l)k3 
~2   +      (s+1) 


(s+1) 


Let  s-*- 1  and  combine  the  non-zero  terms: 

k. 


v  1        2 

K3  "  s(s+l)    (s+1) 


but  k2  =  -1,  hence: 


k3  "  s(s+l)  ' 


Let  s  =  -1,  therefore  the  coefficient,  k3,  is:   k3  =  -1. 
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The  denominator  of  the  unknown  coefficient  will  always 
be  cancelled  when  the  above  procedure  is  used.   The  resulting 
partial  fraction  expansion  is: 

111 


C(s)  = 


1  _ 
s(s+l)2    s    (s+1)2 


1 

Ti+rr 


Example  2-4 


C(s) 


s+2 


(s+1) (s+3) 


TsTiy 


(s+3) 


(s+3) 


(s+3) 


(2-15) 


As  before: 


(s+l)C(s) 


-1+2 


s=-l    (-1+3) 


and 


k2  =  (s+3)JC(s) 


s=-3 


-3+2 
-3+1 


To  find  k3,  multiply  each  side  of  equation  2-15  by  (s+3)  : 


(s+3) ^ (s+2) 
(s+3)3(s+l) 


(s+3)2kx    (s+3)2k2    (s+3)2k3    (s+3)2k4 
(s+1)   +   (s+3)3   +   (s+3)2   +    (s+3) 


Combine  all  terms  that  do  not  go  to  zero  as  s->-3. 

k. 


s+2 


(s+3)  (s+1)    s+3 


but  k~  =  sr,    hence 


v  s+2-s/2-l/2    l/2(s+3) 

K3    (s+1)  (s+3)  "  (s+1) (s+3) 


Therefore : 


k3  ~  2(s+l) 


s=-3 


1 

'4 
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The  same  procedure  is  used  to  find  the  coefficient,  k. ,  i.e., 
(8+3)  (8+2)    (s+3)kl    (s+3)k2    (s+3)k3    (s+3)k4 

+   XT      +   7T      + 


(s+3)3(s+l)      s+1       (s+3)3     (s+3)2      (s+3) 
Let  a-»—  3  and  combine  all  non-zero  terms. 

s+2         k2      k3 

k,  = 


4    (s+3)2(s+l)    (s+3)2    s+3 
But  k2  =  y  and  k^  =  -j,  hence 

It   =  s+2-s/2-l/2+2s/4+3/4  .     (s+3) 
4       (s+1) (s+3)3        4(s+l) (s+3)2 

Let  s  =  -3,  therefore  the  coefficient,  k.  ,  is:   k.  =  —5- 
The  partial  fraction  expansion  is: 

0(.)  =  5±2  .  vs     +  _vi_,  +  ^izi_  +  zizs 

(s+1)  (s+3)3    <s+1>    (s+3)3    (s+3)2    (s+3) 
In  the  preceding  two  examples ,  the  algebraic  manipula- 
tions required  in  the  evaluation  of  the  multiple  root  coeffi- 
cients are  adaptable  to  the  digital  computer.   The  coefficients 
that  were  evaluated  in  the  classical  manner  may  be  determined 
by  associated  polynomials,  specifically  by  the  use  of  Theorem 
1-2  and  corollary  1-3.   The  combined  operation  is  detailed 
below  for  two  possible  cases,  i.e.,  when  the  repeated  roots 
are  of  order  two  and  three. 
Given  the  function: 

m 

b.s1 


i=0  1 N(s) 

(s-P;L)  (s-p2)  .  .  .  (s-pm)  .  .  .  (s-pn_2) 


— — ; 1  =  DiiT     (2_16) 
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The  partial  fraction  expansion  of  C(s)  is 


C(s)  =  ™±K   + 


+ 


(s-Pm) 


m 


2    (s-pm)    (s-p1) 


+  .  .  .  + 


n-2 


(s-p„  0) 
n-z 


(2-17) 


The  coefficients,  k,,  k2 ,  .  ..,  k  _~,  are  found  using  the 
procedure  outlined  in  case  (1)  above.   The  multiple  order 

coefficients,  k  ,,  and  k   are  found  by  the  following  method. 

m+i      m  -1  3 


k™.Li  =  ts-pm)  C(s) 

m+1      rm 


N(s) 


s=pm   D(s)/(s-pm) 


(2-18) 


s=p 


m 


The  denominator  of  equation  2-18  can  be  expressed  using  equa- 
tion 1-19  as: 


D(s)     nr2  h   ?  j-(h+2).  ,   i 

T  =      L  s    I  Pm        A  .  (p  ) 

/     \  *•        un  .  t,«  m        I'm 

(s-p  )     h=0  j=h+2  J 

Ivaluating  at  s  =  p  ,  the  coefficient  k  ,,  equals: 

3        *m  m+1   ■• 

N<Pj 

,  m 


(2-19) 


n 


(2-20) 


m+1   n-2  . 

v  n  v      1~ (h+2) _  ,   x 

Z  P™  L         P«,        A.  pM 

v  L~  ^m  .  t  ,  ->   m        i  ^m 

kh=0  j=h+2  J 

In  order  to  find  the  coefficient,  k  ,  the  algebraic  method 

m        ^ 

shown  in  examples  2-3  and  2-4  must  be  employed,  that  is, 
multiply  each  side  of  equation  2-17  by  the  denominator  of  k 

,     ,    N(s)    (S-Pm)km+1    (s"Pm)km 


m 


m 


s"Pj 


<s-pJ 
^m 


^7 


kl  kn-2 

(S-Pl}  '"    (s-pn-2}J 


(2-21a) 


Let  s->p   and  combine  each  term  that  does  not  go  to  zero, 
m  J 
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m+1 


N(s)-km+1  [ (s-p.^  (s-p2) . . . (s-pn_2>  J 
(s-pm) (s-p1) . . . (s-pn_2) 


(2-21b) 


s=p 


m 


The  numerator  of  equation  2-21b  will  always  contain  a  fac- 
torable term,  (s-p  ) .   Since  (s-p  )  cancels  in  the  combination 

^m  ^m 

of  the  two  terms  which  individually  become  indeterminate  as 
s^-pm,  the  combination  is  finite.   Hence: 

N (s) "km+l (s"pl) (s"p2) ' ' ' (s-pn-2} 


m 


k   = 

m 


(s-p  ) 
■^m 


(s-p1) (s-p2) • • • ^s~Pn_2^ 


(2-22) 


s=p 


m 


Example  2-5 

Expand  by  partial  fractions : 


C(s)  = 


s(s+l) 


N(s) 
DliT 


Ts+iy 


(s+1) 


Solution: 


Let:   p,  =  0;  pw  =  -1. 

ll      •  m 


Form  the  associated  polynomials  of  the  denominator,  D(s): 


A3(u)  =  1 


A2(u)  =  2  +  u 


A,  (u)  =  1  +  2u  +  u' 


Applying  equation  2-7,  the  coefficient  k,  is: 

_    N(0)   . 
Kl  ~  A1(0)+0  "  -1 

Applying  equation  2-20,  the  coefficient,  k 


m+1 


k3  is 


k   -  N(-l) 

K3  "  -A3(-l)  (-l)  +  [A2(-l)  +  (-l)A3(-l)] 


k3  "  -1+1-1  "  _1 
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Applying  equation  2-22,  the  coefficient,  k   =  k2  is 

l-(-l)s 


s+1 


=  -1 


S=  -1 
The  partial  fraction  expansion  of  C(s)  is: 


C(s)  = 


s(s+l)2    s    (s+1) 


1    ,  ,  1 
2  +  (s+1) 


which  agrees  with  example  2-3. 

When  the  expression,  C(s),  contains  a  third-order  pole 
such  that: 


C(s)  = 


m 

I      b.s1 
i=0  x 


N(s) 


(s-pm)  (s-p1) (s-p2) . . . (s-pn_3) 


(2-23) 


The  expanded  form  is: 

-,,      ,   N(s)    _   km+2  km+l               km                  kl                   , 

c(s)  -^-^pT3  7IT77    "^7  +  ^I^ +'"  +^^7 

*-     TYl  ■'-    TT\ 


n-3 


m 


m 


(2-24) 

As  above,  the  coefficients,  k..  ,  k2 , .  .  .  ,  k  _-,  are  found  using 

the  procedure  of  case  (1).   The  multiple  order  coefficients, 

k  xn,  k  ,, ,  and  k   are  found  by  the  following: 
m+2    m+1       m  ■*  3 


km+2 


=  (s-pJJC(s) 


m 


N(S) 


s=pm   D(s)/(s-pm) 


(2-25) 


s=p 


m 


Using  equation  1-21,  the  denominator  of  equation  2-25  can  be 
expressed  as: 


D(s) 


-   I   s^1  I       [j-(h+l)]PJ-(h+2)A,(p  )    (2-26) 
i-  i  ■    t  .  *  m         i      m 


(s-p  r   h=l     j=h+2 


m 
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Evaluating  as  s  =  p  ,  the  coefficient  km+2  is: 
k 


N(Pj 
m 


wm+2        n-2 


I      P^M    j3-(h+l)lpJ-<h+2>A.(Pm, 


(2-2-7) 


h=l      ll1        j=h+2 

The   coefficients,    k    ,,    and  k    ,    are   evaluated   as    in   the 

m+i      m 

second-order  example  above.   These  coefficients  are  found  to 
be : 

N (s)~km+2 (s~Pi)  ^s_p2^  *  * '  ^s-pn-3^ 

(s-p  ) 
^m 


Lm+1        (s-p,) (s-p2) • • •  ^s"pn_3^ 


S  =  D 


m 


(2-28) 


and 


N(s) -k   0 (s-p, ) . . . (s-p„  ,) -k   , (s-p) . . . (s-d_  *) 

m+2    ri       L  n-  J   m+i    m        n-J 


k  = 
m 


(s-p  ) 
-m 


(s-px) (s-p2. . . (s-pn_3) 


S  =  D 


m 


(2-29) 


Example  2-6 

Expand  by  partial  fractions : 


C(s)  = 


s+2 


N(s) 


3    D(s)    s+1 


+ 


(s+3)     (s+3) 


(s+3) 


(s+1)  (s+3) 
Solution : 

Let:   p,  =  -1;  p   =-3 

^1  MTl 


Form  the  associated  polynomials  of  the  denominator,  D(s): 
D(s)  =  s4  +  10s3  +  36s2  +  54s  +  27 
A4(u)  =  1 

A3(u)  =  10  +  u 
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A2(u)  =  36  +  lOu  +  u 

Ax(u)  =  54  +  36u  +  10u2  +  u3 

Applying  equation  2-7,  the  coefficient,  k. ,  is: 

N(-l) 


kn  = 


A  (-!)+(-!) A? (-l)+(-l)2A_(-l)+(-l)3A4(-l) 


where:   N(-l)  =  1; 


Hence 


A,(-l)  =  27;  A2(-l)  =  27;  A  (-1)  =  9;  A,(-l)  =  1 


k,  = 


'1  "  27-27+9-1  "  8 

Applying  equation  2-27,  the  coefficient,  k2 ,  is: 

k   -  N(-3) 

K2    (-3)A4(-3)+[A3(-3)+2(-3)A4(-3)] 


where : 


Hence : 


N(-3)  =  -1;  A3(-3)  =  7;  A4(-3)  =  1. 


k   =     -1     -  4 

K2        -3+7-6    2 
Applying  equation  2-28,  the  coefficient,  k~,  is 

s+2-l/2(s+l) 


k4  = 


k4  =  -1/8 


y       -                     s+3       _  1/2 

3         s+1        s+1 

s  =  -3 

k3  =  -1/4 

equation  2-29,  the  coefficient,  k. ,  is: 

s+2-l/2(s+l)+l/4(s+l) (s+3) 

(s+3)2 

1 

s+1 

4(s+l) 

s  =  -3 
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CHAPTER  III 
Polynomial  Sensitivity  Theory 
3.1:  Introduction 

In  this  chapter  several  simple  algebraic  operations  in- 
volving polynomial  sensitivity  are  investigated.   Specif- 
ically, associated  polynomial  theory  is  used  to  study  the 
relationships  between  variations  of  polynomial  coefficients 
and  polynomial  roots.   Both  large  and  incremental  variations 
are  considered.   As  is  shown  in  the  next  two  chapters,  the 
results  derived  from  this  section  have  very  useful  aoplica- 
tions ,  particularly  in  the  field  of  sensitivity  analysis 
(Chapter  Four)  and  in  self-adaptive  control  systems  (Chapter 
Five).   In  order  to  obtain  these  results,  a  basic  ortho- 
normality  relationship  between  the  set  of  associated  polyno- 
mials and  the  roots  of  the  defining  polynomial  is  defined. 
In  addition  to  being  used  in  the  derivation  of  sensitivity 
expressions,  the  orthonormality  relationship  can  be  used  to 
solve  a  variety. of  linear  problems  that  involve  matrix  in- 
version. 
3.2:  Coefficient  Sensitivity 

The  effect  of  variations  of  a  root  on  the  coefficients 
of  a  polynomial  can  be  expressed  by  the  following  theorem: 
THEOREM  3-1 

The  variations  of  the  coefficients  of  a  polynomial  with 

respect  to  the  variation  of  a  root  may  be  determined  using 

polynomials  by  the  relationship: 

Aa  ._, 

£-J =  -A.(pi)  :  i,j  =  1,2,. ..,n 

pi       J 
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This  result  holds  for  large  variations  when  p.  is  a  simple 

root. 

Proof: 

Let  a  monic  polynomial,  F(s) ,  of  degree  n  with  n  distinct 
roots  be  expressed  as: 

^ /  \     n        n-1  v     i        i 

F(s)=s   +  a  _,s    +...+an=  I      a.sJ;a  =1 
»"*  u    j_0   3       n 

(3-1) 
By  the  fundamental  theorem  of  Algebra,  the  polynomial  may 
also  be  written  as : 


F(s)  =  7X  (s_ph^  =  f (Pl'P2/# * ' ,pi 


, .  .  . ,pn)     (3-2) 


Let  Ap.  denote  a  change  in  the  root,  p..   The  polynomial, 
f(P-wP?/'''/P-+Ap.  ,...,p  )  ,  may  be  expressed  in  a  Taylor  se- 
ries  expansion  about  the  point  (p, ,p2 , . . . ,p. , . . . ,p  )  by: 

f  (Pi  >T?2  ' '  '  *Pi+APi'  *  *  *Pn^  =  f  ^Pl/P2  '  *  *  *  /Pi'  *  '  *  /Pn^  + 

9f (p, ,p~ , . . . ,p.  , .  .  . ,p  ) 

± £ i £L_  Ap.  + 

3pi  *i 

2 

3  f  (Pi  /Po>  •  •  •  /P-!  /  •  •  •  /Pn)         o 

1/2  X      Z    2 1 2-  (APi)2  + 

3Pi 

...  +  i  znf(p1,P2"'"Vi'-'"pn)    tK     ;n    ,,  „ 

E7  — n (Api}      (3"3) 

3p. 

However,  from  the  factored  form  given  in  equation  3-2: 

af(Pl,P2,...,Pi,...,Pn)     F(s) 

^T  =  "  T1RTT       (3"4) 
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and: 

9  f (PirP2' • ' • 'Pi ' • • * 'Pn) 

"~7   n~~ 
3Pi 


=  0  :  n  ?    1        (3-5) 


Equation  3-5  applies  when  p.  is  a  simple  root  only.   Hence 

equation  3-3  becomes: 

f (p1,p2,. . . ,Pi+Api,. •• ,pn)  =  f (p1,p2,. •• ,pn)  -  Api  |i|I 

^i 

(3-6) 

and: 

AF  =  -Ap.  £l£i  (3-7) 

1  s-pi 

where : 

AF  =  f (p1,p2, . . . ,pi+Api/ . . . ,pn)  -  f (p1,p2/ . . . ,pn)  (3-8) 

The  expression  AF  in  equation  3-8  may  be  expressed  in  terms 
of  the  changes  in  the  coefficient,  Aa,  ,  of  the  original 
polynomial,  F(s).   Note  that: 
f (p1,p2, . . . ,pi+Api, . . . ,p  )  ■  (s-p1)  ...  (s-pi+Api)  ...  (s-pn) 

(3-9) 
Substituting  equations  3-2  and  3-9  into  equation  3-8  and  sub- 
tracting the  coefficients  of  equal  Dower  of  s  reveals  that: 

n-1    n  -i-1 

AF  =  AaA  +  Aa,s  +  ...  +  Aa   ,s     =   )  Aa.  ,s:_   (3-10) 

0      1            n-1        ,L,  -i-l 

j«l  J 

Using   Theorem   1-2,    i.e.,; 

^=       I      A. (p. 
'"Pi         i  =  l       ^    Pi 


liSLm        I       XApA)B*-1  d-17) 

s-p.  ,»-       i    FL 


40 


Equation    3-7    is    rewritten   as: 


n 


,3-1 


AF   =    -Ap.       I      A.  (t>.)s3~x  (3-11) 

1    j=l      :      1 


Equating  equations  3-10  and  3-11: 


AF  =  J      Aa.  1sj"1  =  -APi  J      A  (pi)sj'1      (3-12) 


j=l    J  *  -  j-1 

Equating  the  coefficients  of  equation  3-12,  it  follows  that: 

Aa.  x  =  -ApiA.(pi)  :  i,j  =  1,2,. .  .  ,n        (3-13) 

Thus  -A. (p.)  may  be  characterized  as  the  change  in  the 
coefficient  a._,  per  unit  change  in  the  root,  p..   Note  that 
equation  3-13  holds  in  the  large,  i.e.,  for  arbitrary  changes 
in  Ap. ,  when  the  root  is  simple.   Similar  precise  equations 
may  be  derived  for  repeated  roots.   Equation  3-13  is  useful 
for  determining  the  change  in  the  coefficients  of  a  polyno- 
mial when  there  is  a  change  in  a  simple  root.   Equation  3-13 
cannot  be  used  to  determine  the  change  in  root  locations 
caused  by  changes  in  any  coefficient  of  a  polynomial  since 
a  change  in  general  can  be  expected  to  change  the  location 
of  several  roots.   Thus  the  inverse  of  equation  3-13  is  not 
meaningful . 
Corollary  3-1 

Taking  the  limit  of  equation  3-13  as  the  change  in  the 

root,  p.,  goes  to  zero,  then: 

8a  .  , 
a3    =  -A.(p.)  :  i,j  =  1,2,. ..,n  (3-14) 

dPj_        D   ! 
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Example  3-1 

Given  the  polynomial: 

F(s)  =  s3  +  6s2  +  lis  +  6  =  (s+1)  (s+2)  (s+3) 
The  associated  polynomials  are: 

A3(u)  =  1 

A2  (u)  ■  6  +  u 

A1(u)  =  11  +  6u  +  u2 
Let  the  roots  of  F(s)  be  defined  as: 

p1  =  -1;  p2  =  -2;  p3  =  -3 

Assume  p,  is  decreased  by  8  to  a  new  value  p,  +  Ap,  =  -9. 
Hence : 

APl  =  -8 
The  original  polynomial  is  changed  to: 

F1(s)  =  s3  +  14s2  +  51s  +  54  =  (s+9)  (s+2)  (s+3) 
By  inspection,  the  change  in  the  coefficients  are  seen  to  be: 
AaQ  =  48;  Aa,  =  40;  Aa2  =  8;  Aa3  =  0 

In  order  to  verify  theorem  3-1  numerically,  the  changes 
in  the  coefficients  can  also  be  found  by  use  of  equation  3-13 
For  i  =  1 : 

Let  j  =  3,  or: 

Aa2  =  -A3(Pl)APl 

numerically : 

Aa2  =  -(1) (-8)  =  8; 
which  agrees. 

Let  j  =  2,  or: 

Aal  "  "A2(Pl)APl 
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numerically : 


Aa;L  =  -[6+(-l)]  (-8)  =  40; 


which  agrees. 

Let  j  =  1,  or: 

Aa0  =  "Al(pl)Apl 
numerically: 

AaQ  =  -[ll+6(-l)+(-l)2] (-8)  =  48; 
which  agrees. 
3.3:  The  Basic  Orthonormality  Relation 

In  order  to  derive  the  inverse  of  Theorem  3-1,  that  is, 
the  variation  of  a  root  of  a  polynomial  with  the  variation 
of  its  coefficients,  the  orthonormality  relationship  of  asso- 
ciated polynomials  must  first  be  derived.   Prior  to  this 
derivation,  Theorem  3-2  is  introduced. 
THEOREM  3-2 

The  derivative  of  a  polynomial  may  also  be  found  from 
its  associated  polynomials  by: 

SiSL.        m      j    j-1   (u) 

du  j=l       3 

Proof: 

It  was  shown  in  the  proof  of  Theorem  1-2  that: 

F(S)-F(U)  =   f   s^A.U)  (1-14) 

Let  s->u,  then 


n 


,.   F(s)-F(u)  ,.    r  1-1,  /  \          / ->  1c\ 

lim  — ! — J — -  =  lim   )  sJ   A.  (u)          (3-15) 

s-u          . L,  3 

s-*-u  s-*u  j  =  l  J 
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or 


n 


F' (u)  =  I      u?    1A. (u) 

j  =  l      ^ 


(3-16) 


Example  3-2 


Given  the  polynomial: 
F(s)  =  s4  +  10s3  +  35s2  +  50s  +  24  =  (s+1)  (s+2 )  (s+3)  (s+4) 
The  associated  polynomials  are: 
A4(u)  =  1 

A3(u)  =  10  +  u 

A2(u)  =  35  +  lOu  +  u2 

A1(u)  =  50  +  35u  +  10u2  +  u3 

The  derivative  of  the  polynomial,  F(s),  is  expressable  in 

two  ways : 

(1) :   Classical  Method: 

dFis)  =  4s3  +  30s2  +  70s  +  50 
ds 

(2):   Associated  polynomial  method: 

Applying  Theorem  3-2; 

d||^-  =    [A1(u)+A2(u)u+A3(u)u2+A4(u)u3]u=s 

=  50  +  35u  +  10u2  +  u3 

+  35u  +  10u2  +  u3 

2     3 
+  lOu   +  U 

+  u3 


d^s)  =  50  +  70u  +  30u2  +  4u3 


u=s 
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THEOREM  3-3 

The  basic  orthonormality  relationship  between  the  set  of 
associated  polynomials  and  the  roots  of  the  defining  polyno- 
mial is: 

n  Ai(p  ^p1"1 

iil   F"(Pi'       kjj,k=l,2 n 

where  6,  .  is  the  Kronecker  delta  function  which  vanishes  ex- 
cept when  the  indices  are  equal. 
Proof: 

The  derivative  of  a  polynomial  is  expressable  using 

Theorem  3-2  as: 

n    .  , 
F'  (u)  =  I      u^'-V  (u)  (3-16) 

j-1      : 

If  p.  and  p,  are  roots  of  the  defining  polynomial,  F(s), 

then  it  can  be  shown  that,  if  p.  is  a  simple  root  of  F(s) : 

n 


0  =  I      p  riA  (p.)  (3-17) 

i=l   K     3      1 


J- 

From  equations  3-16  and  3-17,  it  follows  that: 
n 


I   P^TVpi) 


j  pi      0  when  i  j*   k 

^ =  <  ,   .    •    .  (3-18) 

n     .  ,  1  1  when  i  =  k 

.1   P^^A  (p.) 

3=1 


Equation  3-18  may  be  rewritten 

j-l 


n  A.  (p.)  p. 

V    3  )      * -6..  (3-19) 

S    F'  (p. T       ik   ,    ,  ~ 

=1      *i  i,k  =  1,2 , . . . ,n 


where  6.,  is  the  Kronecker  delta  function, 
lk 
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Let  [a]    and  [3]  denote  the  matrices  who  elements  are: 

A  (p.) 

a  .   .  =  — 4—, r  j  =  1,2,. ...n;  i  =  1,2,.  .  .  ,  n     (3- 20a) 

ID    F  (Pj_) 

and 

3  -k  =  p,-1      j  =  1,2,. ..,n;  k  =  1,2,.  .  .  ,n     (3-20b) 

If  I  denotes  the  identity  matrix,  then  equation  3-19  may 

be  written: 

n 

[a]  [6]  =  I      a.  .3.,  =  6.,  =  I  (3-21) 

=1   J-D  j k     ik 

Since  3 -k  is  a  Vandermonde  matrix,  it  is  non-singular  and: 

n 
[3] [a]  =  I      a    .a. .  =  I  (3-22) 

iil   kl  ^ 

Hence : 

n   A. (p.) p. 

j    J   )   ^ =  6  (3-23) 

i=l   F  ipi;       K:   j,k  =  1,2,. .  .  ,n 

Theorem  3-3  may  also  be  verified  by  the  following  example: 

Example  3-3 

F(s)  =  s3  +  6s2  +  lis  +  6  =  (s+1)  (s+2)  (s+3) 

Let: 

Pi  =  -1'    P2  =  "2;  P3  =  ~3, 
The  associated  polynomials  of  the  given  polynomial  are: 

A3(u)  =  1 

A2(u)  =  6  +  u 

9 
A, (u)  =  11  +  6u  +  u 

The  derivatives  of  the  polynomial  evaluated  at  the  three 
roots  are  found  using  equation  3-16  to  be: 
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F'(p.)  =  Al(Pi 


P'  (p±)  =6+5 


F'  (p2)  =  3  +  4 
F' (p3)  =2+3 


Forming  the 
A. 


[a..]  = 


F^ 
A, 


F" 
A, 


a]  and 
Pl)   A2 


Pi)   F' 

pTT  r 


P3)   A2 


pTT  f^ 


3] 


Pi 
p2 


p2 
P3 


+  A2(Pi)p.  +  A3(Pi)pi 
-1)  +  1  =  2 
-2)  +  4  =  -1 
-3)  +  9  =  2 
matrices  according  to  equation  3-20 


A3(p1) 


F' (Pl) 
A3(p2) 


A3(p3) 


f-T^Tj 


6/2       5/2       1/2 
3/-1    4/-1    1/-1 
2/2      3/2      1/2 


Or 


[a.  .]    = 
ID 


3      5/2      1/2 
-3         -4      -1 
_1      3/2      1/2  J 


i<V  = 


0 


P 


0 


111 

Pl     P2     P3 


111 

-1      -2      -3 
14         9 


It   can  be    shown   that 


Likewise : 


[a±jl[pjk]  =  I 


[3.,  ]  [a.  .]  =  I 
jk   13 


(3-21) 


(3-22) 
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Expanding   the    terms    in   equation    3-22: 


0         0 


P- 


111 
Pi      P2      p3 


Lpl 


'3j 


Hence : 


A1(P3 
FTTpT 


0    Al(pl}  0    Al 

Pl      F'  (Pl)  P2      F1" 


Al(Pl 


A]_(p2 


F« 
A„ 


F 
A^ 


F' 


PX)  A3 

pTt  r 

p2)  A3 

pT"  f- 

p.)  A^ 


pT7    f1 


Po) 


0   Al 

+  p3    f 


P? 


P^T 


p2) 


P^" 

p.) 


p?j 


=  1 


p.) 


pV> 


=  1 


Or  in  matrix  notation, 


likewise : 


and: 


[sj=i,k1[ai,j=i!  -  L 


[gj=2,k][ai,j=21  "  * 


In  addition,  the  following  relationships  also  hold: 

[ej  =  l,k]K,j  =  2]  =  ° 


(3-24a) 

(3-24b) 

(3-24c) 

(3-24d) 

(3-25a) 
(3-25b) 


[ej  =  3,k]^ifj  =  2]  =  0 


(3-25c) 


The  above  matrix  products,  equations  3-24  and  3-25  can 


be  expressed  explicitly  as  stated  in  Theorem  3-3,  i.e., 

n   A.(p.)p.k_1 

I         3r, )_    * =  6VJ  (3-23) 


i^-np? 


kj 


48 


The  orthonormality  relation  given  by  equation  3-23  is  a 
basic  result  which  permits  one  to  solve  a  variety  of  linear 
problems.   As  has  been  shown  in  references  1  and  11,  this 
orthonormality  relation  permits  the  theoretical  investiga- 
tion of  inverse  matrix  problems,  such  as  finding  the  inverse 
to  the  Vandermonde  Matrix  in  a  form  suited  for  numerical  cal- 
culations . 
3.4:  Root  Sensitivity 

The  effect  of  variations  of  the  coefficients  of  a  polyno- 
mial on  any  particular  root  is  expressable  by  the  following 
theorem: 
THEOREM  3-4 

The  variation  of  the  roots  of  a  polynomial  with  respect 
to  the  variations  of  a  coefficient  of  the  defining  polyno- 
mial is  determined  by  the  relationship: 

9pi     pi:       j  =  0,1,.. . ,n-l 

3T7  =  "FHTp.)  :  .   ,  9 

j  "i     l  =  1 , 2  , .  .  .  ,n 

Proof: 
Corollary  3-1  established  that: 

da  ._, 

a-l~  =  -A.  (p.)  :  i,j  =  1,2,.  .  .  ,n 

dp^      j      1 

In  addition  the  orthonormality  relation  of  Theorem  3-3  is: 

n   TV  /   \   k-1 

I      Ai(Pi)Pi     =  «ki  (3-23) 

i=l    F' (p.)         K: 

Let  [y]    and  [r\]    denote  the  matrices  whose  elements  are: 

y..    =   A. (p.)   i  =  l,2,...,n;  j  =  1,2,. . . ,n     (3- 26a) 
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and  Jc.1 

pi 
nik  =  F' (p. )   i  =  l»«,...,nj  k  =  1,2,. ..,n     (3-26b) 

Using  the  matrix  notation  introduced  above,  equation  3-23 
can  be  rewritten: 

[y] [n]  =  I  (3-27) 

Hence,  the  matrix,  y,  (whose  elements  are  A. (p.))  is  the  in- 

k-1 
verse  of  the  matrix,  n,  (whose  elements  are  p.    /F'(p.)) 

given  i,j,k  =  l,2,...,n.   In  a  similar  manner,  it  has  been 

shown    that  the  matrix  with  elements  9p./9a.  ,  is  the  in- 

verse  of  the  matrix  whose  elements  are  8a._,/8p.,  that  is, 

if: 

and 


[A]  =   v-i  (3-28a) 

da . 


8ak-l 
[B]  =    .K  X  (3-28b) 

8p. 

Then: 

[A] [B]  =  I  (3-28c) 

but  from  corollary  3-1: 

[B]  =  -[y]  (3-28d) 

Hence,  from  equations  3-27  and  3-28d: 

[A]  =  -[n]  (3-28e) 

Writing  out  the  terms  of  the  matrices,  as  given  by  equations 
3-28a  and  3-26b  yields: 
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dp±  dp1 


9a   9a, 
o    1 


9p. 


9a, 


dp 


n 


L8a0 


9p. 


9a 


n-1 


9p 


n 


9a 


n-]J 


P- 


F'  (Pl)  F'  (P]_)  * 


P- 


FT5J) 


n 


F'  (p  ) 


n-1 


F' (Pl) 


n-1 
_Pn 

F'(p  ) 


(3-29) 


Thus  in  general : 


9pi     p^  i  =  1,2,. . .  ,n 

9a.  =  ~F' (p. )       L  ,        , 
D         i    D  =  0,1,... ,n-l 


(3-30) 


-th 


Corollary  3-2 

The  total  change  in  the  value  of  any  root  of  an  n 
order  polynomial  due  to  an  incremental  change  in  any  or  all 
coefficients  of  the  polynomial  is  found  by: 


dp.  = 


9p .        9p .  9p 


9a 


0 


where : 


da0  +  9a. 


9p, 


9a 
m 


da,  +  . .  .  + 


m 
pi 

T  (pi) 


da 


9a   ,    n-1 
n-1 


(3-30) 


Substituting  and  factoring: 


dp.  = 


dan  +  p . da,  +  ...  +  p 


n-1 


da 


n-1 


-F' (p±) 


(3-31) 


3.5:  Root  Solving  Using  Associated  Polynomials 

Equation  3-31  of  Corollary  3-2  may  be  used  in  a  root 
solving  process  if  a  polynomial  with  known  simple  roots  can 
be  found  near  the  unknown  polynomial.   Since  equation  3-31 
deals  with  differentials,  the  process  is  limited  to 
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polynomials  whose  roots  are  close  in  value  to  the  known  polyno- 
mial, i.e.,  small  variations  in  the  coefficients  between  the 
two.   The  process  is  summarized  as: 

(1) :   Find  a  polynomial  of  equal  order  with  known,  non- 
repeated,  roots,  such  that  the  difference  in  the  coefficient 
values  are  "small".   The  accuracy  of  this  root  solving  proc- 
ess is  determined  by  the  value  of  the  variations  in  the  coeffi- 
cients . 

(2) :  Using  the  roots  of  the  known  polynomial  and  the 
values  of  the  change  between  the  coefficients  of  the  known 
and  unknown  polynomials,  apply  equation  3-31. 

(3):   The  change  in  the  known  roots  determines  the 
values  of  the  roots  of  the  unknown  polynomial. 

The  following  examples  demonstrate  how  this  root  solv- 
ing process  may  be  employed. 
Example  3-4 

The  roots  of  a  quadratic  polynomial  are  easily  deter- 
mined by  the  quadratic  formula.   However,  such  analytical 
expressions,  in  general,  do  not  exist  for  higher  order  poly- 
nomials.  This  example  will  use  a  quadratic  polynomial  as 
an  illustration  of  the  root  solving  method  available  with 
associated  polynomials. 
Given  the  quadratic  polynomial: 

G(s)  =  s   +  4.7s  +  5.3 
In  order  to  find  the  roots  of  G(s) ,  note  that  the  polyno- 
mial ,  F (s)  ,  is  : 

F(s)  =  s2  +  5s  +  6  =  (s+2)  (s+3) 
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The  variations  in  the  coefficients  of  F(s)  from  the  qiven 
polynomial,  G(s) ,  are  small  enough  such  that  the  process  may 
be  employed  with  negligible  errors. 

Let  the  roots  of  F(s)  be  defined  as: 

pl  =  ~2;  p2  =  "3  * 
Evaluating  the  denominator  of  equation  3-31: 

F'  (s)  =  2s  +  5 
Therefore : 

-Ff  (Pl)  =  -1;  -F' (p2)  =  +1 

The  variations  of  the  coefficients  are  identified  as: 
daQ  =  -.7;  da.  =  -.3;  da2  =  0. 

The  change  in  the  root,  p.,  of  F(s)  is: 

,       .7+(-2) (-.3)     , 
dp1  =  zi •1 

Hence:   p,  of  G(s)  is  -2  +  .1  =  -1.9. 

The  change  in  the  root,  p2  ,  of  F(s)  is: 

x*           -.7+(-3)  (-.3)     2 
dp2  =  n .  2 

Hence:   p2  of  G(s)  is  -3  +  .2  =  -2.8. 

The  factors  of  the  polynomial,  G(s),  are: 

(s+1.9)  (s+2.8) 
As  a  check  on  the  above  results : 

(s+1.9)  (s+2.8)  =  s2  +  4.7s  +  5.32 
Example  3-5 
Given  the  third  order  polynomial: 

G(s)  =  s3  +  5.85s2  +  10.44s  +  5.39 
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To  find  the  roots  of  G(s),  note  that  the  polynomial,  F(s), 
is  : 

F(s)  =  s3  +  6s2  +  lis  +  6  =  (s+1)  (s+2)  (s+3) 
Let  the  roots  of  F(s)  be  defined  as: 

pl  =  ~1?  p2  =  ~2;  p3  =  "3* 
The  denominator  of  equation  3-31  is: 

F*  (s)  =  3s2  +  12s  +  11 
Therefore : 

P1 (p£)  ■  3  -  12  +  11  =  2 

F'  (p2)  =  12  -  24  +  11  =  -1 
F1 (p3)  =  27  -  36  +  11  =  +2 
The  variations  of  the  coefficients,  [a.  of  G(s)  -  a.  of  F(s)], 


are : 


da~  =  -.61;  da,  =  -.56;  da-  =  -.15 


The  change  in  the  root,  p.,  of  F(s)  is: 

da0+p1da1+p12da2 
dpl  "  =F7p^l 


,              (-.61)+(-l) (-.56)  +  (-lP  (-.15)     . 
dpl =2 "  -1 

Hence:   p,  of  G(s)  is  -1  +  .1  =  -.9. 

The  change  in  the  root,  p2 ,  of  F(s)  is: 

2 

da,+p2da,+p2  da2 

dp2  "  ^F7!^) 


A«  (-.61)+(-2)  (-.56)-K-2)2(-.15)  _    nq 

ap2  _  -.uy 
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Hence:   p2  of  G(s)  is  -2  -  .09  =  -2.09. 

The  change  in  the  root,  p3,  of  F1 (s)  is: 

da0+p3da1+p32da2 
dp3  =  -F'(p3) 

dp   =  (--6D  +  (-3)  (-.56)+(-3)2(-.15)  =  ^  14 

Hence:   p3  of  G(s)  is  -3  +  .14  =  -2.86. 
The  factors  of  the  polynomial,  G(s),  are: 

G(s)  =  (s+.9)  (s+2.09)  (s+2.86) 
As  a  check  on  the  above  results: 

(s+. 9)  (s+2.09)  (s+2.86)  =  s3  +  5.85s2  +  10.43s  +  5.38 
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CHAPTER  IV 
Sensitivity  Coefficients 

4.1:  Introduction 

In  the  analysis  of  a  system  whose  motion  can  be  described 
by  a  set  of  differential  equations,-  in  addition  to  obtaining 
solutions  as  functions  of  independent  variables,  the  engineer 
would  like  a  knowledge  of  the  variations  of  the  solutions 
with  respect  to  the  parameters  of  the  system.   Since  it  is 
impossible  to  realize  a  physical  system  that  is  identical  to 
a  mathematical  model,  the  influence  which  parameter  varia- 
tions in  the  physical  system  have  on  the  system's  behavior 
is  an  important  condition  for  reliable  design. 

Variations  in  a  circuit  element  will  cause  a  displace- 
ment of  the  poles  and  zeros  of  driving-point  and  transfer 
functions.   As  a  direct  consequence,  these  variations  and 
resulting  displacements  Effect  the  network  response.   One 
method  of  approaching  this  problem  is  to  plot  the  displaced 

poles  and  zeros  and  then  determine  the  new  response.   As  this 

(12) 
method  involves  considerable  work,  it  has  been  shown     that 

the  work  involved  may  be  reduced  by  testing  the  driving-point 

or  transfer  function  to  determine  incremental  changes  in  a 

particular  parameter.   This  method  is  called  "Sensitivity 

Analysis" . 

(13) 
Tomovic  introduced     sensitivity  coefficients  for 

error  and  sensitivity  analysis  of  dynamic  systems.   He  has 
shown  that  determination  of  the  influence  of  parameter  tol- 
erances on  the  behavior  is  based  on  a  knowledge  of  these 
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coefficients.   They  show  the  effects  of  incremental  parameter 
variations  on  the  system's  behavior  by  giving  the  variation 
of  the  system  response  for  any  given  input. 

Systems  which  contain  parameters  that  undergo  large 
variations  are  studied  in  the  next  chapter.   This  chapter  ex- 
presses Tomovic'  sensitivity  coefficients  analytically  in 
terms  of  associated  polynomials.   Tomovic'  work  demonstrated 
how  these  coefficients  may  be  obtained  from  a  computer  model 
and  analytic  expressions  for  them  have  not  been  available. 
Applications  of  sensitivity  coefficients  for  error  and  stabil- 
ity analysis  has  been  amply  covered  in  reference  13  and  is 
not  covered  here. 
4-2:  Theory  of  Sensitivity  Coefficients 

Any  dynamic  system  may  be  represented  by  an  input,  an 
output  and  a  transfer  function  relating  the  two  as  shown  in 
figure  4-1. 


Input  y (t) 


System 


Output  x(t) 


Figure  4-1:   Dynamic  System  Basic  Block  Diagram 

If  the  Laplace  transformation  of  the  input  signal,  Y(s), 
is  equal  to  one,  then  the  response  of  the  system  is  dependent 
on  n  parameters  and  is  defined  by: 

bQ+b1s+. . .+b  s 

X(s,aQ,a1, .  .  .  ,an)  =  -         (4-1) 

aA+a, s+. . .+a  s 
0   1       n 


57 


Tomovic  introduced  the  complex  sensitivity  coefficient, 
V.  (s)  ,  which  can  be  defined  with  respect  to  the  relative 

variation  of  the  parameter  a. ,  as: 

dX(s,an,a, , .  . .  /  a  ) 

V.  (s)  =  "  / 2-  (4-2) 

1  din  a. 

It  follows  from  the  Laplace  integral's  property  of 

differentiability  with  respect  to  the  parameters  that  the 

real  sensitivity  coefficient,  v. (t) ,  can  be  obtained  by  the 

inverse  transformation: 

dX(t,an,a-,  ,.  .  .  , a  ) 
vi(t>  ■  3TH-Si ~  (4-3' 

1 

In  his  text,  Tomovic  continues  from  this  definition  to 
a  method  of  determining  the  sensitivity  coefficients  by 
modeling  on  an  analog  computer.   This  section  uses  asso- 
ciated polynomial  theory  for  the  determination  of  the  sensi- 
tivity coefficients  in  terms  of  the  coefficients  of  the 
system's  transfer  function  which  therefore  makes  them  amen- 
able to  digital  computer  computation. 

The  Laplace  transformation  of  the  response,  equation  4-1, 

may  be  written  as : 

n  n 

I     "»4«    .J.   I     b'«3 

X(s,a.)  ,  J  =  0   3    =   a"   J7°  \ .___ ,     (4-4) 

i     n     .    (s+p, ) (s+p0) . . . (s+d  ) 
r      1       c  l  ^z  n 

I    ais 

i=0   1 
where  the  denominator  has  been  factored.   Expanding  equation 
4-4  into  its  partial  fractions: 

k        k  k 

X(s,a.  )  =  >    J-    y.    +    ,    J"    s  +  ..  .  +  ;   "  .      (4-5) 
l     (s+p1)    (s+p2)  (s +pn) 
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From  equation  4-5,  two  operations  may  be  noted: 

k. 


(1)  : 


(2)  : 


3X(s,a. ) 
i 

3Pi 


3X(s,a. ) 

3k. 


(s+Pi) 


TsTp~T 


(4-6) 


(4-7) 


The  response  is  a  function  of  the  parameters  a.  and  the 

residues ,  k . . 
l 

The  complex  sensitivity  coefficient,  equation  4-2,  may 

be  written  as : 

3X(s,a.  ) 


Vs)  -  ai^r 


n 


=  ai  J 
j-l 


3X(s,a.)  3p.    3X(s,a.)  3k. 


3Pj 


3a. 

l 


3k 


3a 


(4-8a) 


(4-8b) 


3  i 

From  associated  polynomial  theory,  the  variation  of  a 

particular  root  with  the  variation  of  a  coefficient  is  ex- 

pressable  as: 


3£i 

3a. 
3 


J_  mm       -L  f  £   f   •  •  •  f  1 1 


F  (p.)   .  _  n  , 

ri    3  =  0,1, . . . ,n 


(4-9) 


For  notational  purposes,  let: 


3a~ 

l 


3p. 


3p 


=  Q,5 


n 


1'  3a,  "  Q2?  *  *  * ;  3a. 


=  O 


'n 


(4-10) 


l  l 

Using  equations  4-10,  4-6  and  4-7,  the  complex  sensi- 
tivity coefficients,  equation  4-8  may  be  simplified  to: 


n 


V  (s)  =  a   I 
j-1 


-k. 
3 


3k 


Q.  + 


2. 


i  \2   j    (s+p.)  3a. 

(s+p.)    J       *V    1 


(4-11) 
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Equation  4-11  may  be  further  simplified  by  noting: 

9k.    3k.  3p. 

3i1=  3?1  TE1   '  *  -  i'2 n  (4"12a) 

or  using  equation  4-10 , 

3k.    3k. 

TI1   =  3^"  Qi  ;  j  =  1'2'""n  (4-12b) 

i    pj   J 

The  residue,  k,,  in  equation  4-12b  may  be  expressed  in  terms 

of  the  original  equation  for  the  Laplace  transformation  of 

the  response.   From  that  equation,  for  notational  purposes 

let: 

n 
N(s)  =  I   b.s11  (4-13a) 

j-0  : 

n 
D(s)  =  I      a.s3  (4-13b) 

j-0   3 

Form  the  set  of  associated  polynomials  of  the  denominator 

D(s).   From  associated  polynomial  theory,  the  residue  k.  is 

expressable  as: 

N(p  ) 
k   =  — J j  =  1,2,. ..,n       (4-14) 

J       V    /    xm-l_   ,    . 

)   (p.)    A_(p.  ) 
m=l   J        J 

so  that: 

3k        I         N(p  ) 

5-1  =  ^-J   3 V  (4-15) 

3p .    3p . ^   n         , 

■^    ^   y  (p.j^-l  (  } 

m=l   3     n  - 3  y 

Using  equations  4-12  and  4-15,  the  complex  sensitivity  coeffi- 
cient equation  4-11  may  be  written  as: 
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Vi(s) 


a. 

1 


/  "kl   +   3 


N(P1) 


(s+pl)   S+Pl    8pl/   ?    m-1   , 

m=l 


Q2   /  -k2 


(s+P9)  I  s+p0  +  Zp0\     n 


N(p2) 


2     -  2  >  V    m— 1_  , 


m=l 


-   |    n   ,   a   f       *<Pn> 
Ti+PTTI  s+pn    9pni  n      , 

(  I   Pn   W 

v  m=l 


2. 

"n' 


+  .  .  . 


(4-16a) 


Or 


V,  (s)  =  a,   v     l 


*,  (s+ 


pT> 


j-i  ,*T  L  "  *j 


-k 


s+p .    3p.    n 


N(Pi) 


m 


)  p  .    A  (p  . ) 


(4-16b) 


This  expression  is  cumbersome  but  does  give  an  analytic  ex- 
pression for  the  transform  of  a  sensitivity  function. 
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CHAPTER  V 
A  Two  Parameter  Adaptive  System  Using  Associated  Polynomials 

5.1:  Introduction 

A  feedback  control  system  can  be  represented  by  the  basic 
diagram  as  shown  in  figure  5-1,  or  by  some  modification  or 
extension  to  it.   The  plant,  control  and  feedback  systems  are 
composed  of  real  physical  components.   For  purposes  of  anal- 
ysis and  design  the  physical  system  is  usually  modeled  math- 
ematically in  order  that  the  Laplace  Transforms  may  be  em- 
ployed. 


Input 
R(s) 


Control 


Plant 


Feedback 
Controller 


Output 


C(s) 


Figure  5-1.   Feedback  Control  System 

The  characteristic  equation  of  the  system  is  expressable  as: 

1  +  HG(s)  =  0 
where : 

H(s)  is  the  feedback  transfer  function 

G(s)  is  the  direct  or  forward  transfer  function  com- 
posed of  the  individual  transfer  functions  of  both  the  control 
system  and  the  plant. 

HG(s)  is  the  open-loop  transfer  function. 

The  characteristic  equation  determines  the  general  form 
of  the  system's  natural  response  and  its  roots  give  the  time 
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constants  of  the  transient  exponential  decay  factors.   Stab- 
ility of  the  overall  system  is  found  by  factoring  the  char- 
acteristic equation.   The  characteristic  equation  has  real 
coefficients  which  are  algebraic  combinations  of  physical 
parameters,  which  are  real  numbers. 

Basic  feedback  control  theory  usually  assumes  time- 
invarient  elements  for  purposes  of  design  and  analysis.   How- 
ever parameters  of  the  system  may  be  subject  to  variations 
due  to  a  changing  environment  or  similar  phenomena.   If  there 
is  a  constant  loop  gain  or  a  fixed  compensation  scheme  and 
a  system  parameter  were  to  vary,  the  closed  Iood  performance 
could  be  seriously  deteriorated  and  even  system  instability 
occur.   Examples  of  such  a  changing  parameter  may  be  found 
in  different  operations  such  as  the  change  in  the  aero- 
dynamic parameters  of  an  aircraft  as  it  changes  altitudes 
or  the  effect  of  gravity  on  system  parameters  in  the  Ship- 
board Inertial  Navigation  System,  (SINS) ,  of  a  deep  diving 
submarine . 

If  parameter  variations  are  large  and  slow  compared 
with  the  normal  system's  response  times,  it  is  possible  to 
design  for  the  capability  of  continuously  measuring  the 
variations  and  then  changing  compensation  so  that  the  sys- 
tem's performance  criteria  are  always  satisfied.   Such  sys- 
tems which  automatically  exhibit  the  characteristics  of 
adapting  to  variations  in  system  parameters,  are  called 
Self-Adaptive  control  systems. 
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INPUT 


Q* 


COMPENSATOR 
G1(s) 


PLANT 


K 


G2(s)  "  T^W) 


FEEDBACK 
H(s) 


OUTPUT 


Figure  5-2.   Control  System  with  Variable  Parameter  3 


The  parameter  3  of  the  plant  shown  in  figure  5-2,  is 
known  to  vary  slowly  with  time  or  environment.*   This  param- 
eter appears  algebraically  in  the  coefficients  of  the  sys- 
tem's characteristic  equation.   Since  the  dominant  root  val- 
ues are  determined  by  the  coefficients  of  the  defining  equa- 
tion, the  system's  dominant  pole  location  is  a  function  of 
the  value  of  the  variable  parameter  3.   As  the  value  of  3 
varies  so  also  will  the  location  of  the  dominant  root  or 
roots  of  the  characteristic  equation  hence  the  system's  re- 
sponse to  an  input  and  decay  factors  will  change.   If  the 
new  value  of  3  M^y/e^  the  location  of  the  dominant  root  or 
roots  into  the  right-hand  s-plane,  the  system  becomes  un- 
stable . 

A  block  diagram  of  an  adaptive  control  system  is  shown 
in  figure  5-3.   As  before  the  plant  has  a  parameter  3  which 
is  known  to  vary  with  time  or  environment.   The  "identifica- 
tion and  adapter"  block  continuously  measures  the  input  m(t), 


*In  the  analysis  which  follows,  the  assumption  of  quasi- 
stationary  responses  is  made  in  which  parameters  are  assumed 
to  be  varying  slowly  such  that  the  response  is  approximately 
given  by  the  system's  response  with  the  fixed  parameters  re- 
placed by  a  time  varying  function . (23) 
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COMPENSATOR 
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m(t),    ^ 

PLANT 
G2(s)  =  f(3) 
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> 

FEEDBACK 

f 

H(s) 

Figure  5-3.   Self-Adaptive  Control  System 

and  the  output  of  the  plant  in  order  to  identify  the  param- 
eter 3.   When  a  change  in  3  has  been  detected  the  "identifica- 
tion and  adapter"  block  adjusts  the  variable  comDensation 
block  parameter  a  such  that  the  system  is  returned  to  its 
original  specification.   The  design  of  the  "identification 
and  adapter"  block  is  the  critical  portion  of  any  self- 
adaptive  system  and  involves  two  requirements: 

(1)  Some  method  must  be  devised  for  identifying  the 
change  within  the  system's  plant, 

(2)  A  means  of  adjustment  must  be  provided  which  is 
capable  of  counteracting  the  effects  of  the  parameters  varia- 
tions . 

Identification  of  any  change  in  the  parameters  of  a 
plant  have  been  covered  in  references  14  and  15.   References 
16  and  17  are  concerned  with  methods  for  providing  adjustment 
schemes.   This  chapter  shows  how  associated  polvnomials 
theory  may  be  used  in  the  development  of  a  "compensating 
ratio"  between  two  variables  parameters  of  a  feedback  con- 
trol system.   It  shows  that  a  change  in  a  variable  parameter 
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in  the  plant  may  be  compensated  for  by  a  change  in  another 
variable  parameter  such  that  dominant  root  performance  or  be- 
havior is  controlled,, 
5.2:  Basic  Theory 

The  transfer  function  of  a  simple  feedback  control  sys- 
tem as  shown  in  figure  5-1  is  written  as: 

T(s)  - 1  (s>  ■  iwrir  (5"1> 

The  denominator  of  equation  5-1  when  equated  to  zero  is  the 
characteristic  equation  of  the  system.   This  equation  is  a 
polynomial  of  degree  n  and  may  be  represented  as: 

W(s)  =  a  sn  +  a   tS11"1  +  ...  +  an  =  0        (5-2) 
n      n-1  0 

By  the  fundamental  theorem  of  algebra,  the  characteristic 

equation  has  n  roots.   Hence,  equation  5-2  may  be  written: 

W(s)  =  (s+p^  (s+p2)  .  .  ,  (s+pn)  =  0  (5-3) 

Consider  the  coefficients  of  the  characteristic  equa- 
tion, i.e.,  aQ ,  a,,  ...,  a  .   Assume  any  two  variable  param- 
eters a  and  3  appear  linearly  in  any  number  of  coefficients, 
i.e.  , 

a.  =b.a+c.3+d.  ,c  ^ 

3  3  3  3  (5-4) 

:  j  =  0,1,.  .  .  ,n 

where  b,,  c.,  and  d.  are  real  constants.   Note  that  the 
j    3  3 

coefficient  a .  is  a  linear  function  of  a  and  3.   This  repre- 
sentation of  variable  parameters  in  the  coefficients  of  the 

(18) 
characteristic  equation  was  first  introduced  by  D.  D.  Siljak 

(19) 
and  later  expanded  by  G.  J.  Thaler,  et  al .     ,  in  the 

development  of  parameter  plane  theory. 
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From  general  polynomial  theory,  the  change  in  any  root  is 
determined  by  the  change  in  the  coefficients  of  the  defining 
polynomial  and  vice-versa.   The  roots  of  the  characteristic 
equation  can  be  considered  as  functions  of  the  characteristic 
equation's  coefficients.   Hence  the  differential  p.  (a~ ,a,  , .  .  .  ,a  ) 
can  be  expressed  as : 


dPi  = 


3p. 
3i7da0  + 


9p. 

3^*1  + 


+ 


n 


(5-5) 


From  equation  5-4,  the  coefficients,  a.:  j  =  0,1,..., n, 
are  themselves  functions  of  the  variable  parameters  a  and  3. 
Hence  the  differential  of  any  coefficient  is: 


3a. 

da .  =  •*— 
j         3a 


ida  + 


3a. 


33 


IdS 


(5-6) 


Substituting  equation  5-6  into  equation  5-5: 

3a, 


dp 


i  =  3^[3F^a  +  W-**\   +    3^  {TZ~*a  + 


3a. 


36 


<i3 


3p.   3a 


n 


3a 
n 


3a 


da  + 


3a 


(5-7) 


From  equation  5-4 : 

3a 


rl=b, 
3a      j 


3a. 


( j  =  0,1,. . . ,n) 


(5-8) 


33 


1=  c. 


Substituting  equation  5-8  into  equation  5-7  and  then  factor- 
ing da  and  d3 ; 
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I 


3p. 


3p. 


0  8a, 


+  c 


1  8a 


n  9a 


3p. 

i  i 

+. . .+  c 


'n  8  a. 


da  + 


dg 


(5-9) 


"1        "   nj 
From  associated  polynomial  theory,  the  variation  of  a  polyno- 
mial's roots  with  the  variation  of  its  coefficients  is  exnres- 

sable  as : 

J 


da, 
3 


F'tp.) 


j  =  0,1,.  .  .  ,n 


(5-10) 


Using  the  above  equation,  the  differential  of  any  root,  p. , 


may  be  written  as: 


n 


bA+b. p»+. . .+b  p 
-,        0   l-i      n*i 
dp  .  =   „,  >^ — < da 


-F' (P .) 


co+ciV--Vi 

-F'(Pi) 


n 


de 


(5-11) 


In  this  chapter  the  term  "compensating  ratio"  is  used 
for  the  mathematical  relationship  between  two  variable  param- 
eters such  that  the  system's  dominant  performance  or  behavior 
is  controllable.   Since  the  dominant  root  of  a  feedback  con- 
trol system  may  be  either  real  or  a  complex  conjugate  Dair, 
four  different  cases  of  compensating  adaption  are  possible. 
These  cases  are: 

(1)  The  dominant  root  is  real  and  is  to  be  held  con- 
stant, 

(2)  The  dominant  roots  are  a  complex  conjugate  pair 
and  are  to  be  held  constant, 
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(3)  The  dominant  roots  are  a  complex  conjugate  Dair  and 
the  damping  coefficient  is  to  be  held  constant, 

(4)  The  dominant  roots  are  a  complex  conjugate  Dair  and 
the  damped  resonant  or  natural  frequency  is  to  be  held  con- 
stant. 

5.3:  Case  (1) :   The  Dominant  Root  is  Real  and  is  to  be  Held 
Constant 

When  the  dominant  root,  p.,  is  real  and  is  to  be  held 
constant  for  changes  in  the  parameters  a  and  6,  then  the 
differential  of  the  root  must  be  equal  to  zero,  i.e., 

dPi  =  0. 

With  dp.  =  0,  equation  5-11  reduces  to: 


,      C-+C, p. +c0p.  +...+C  p. 
da      0   l^i   2*1 n  l 

b^+b-p.+b^p.  +...+b  p. 
0   l-i   2^i       n^i 


(5-12) 


Where  da/d3  is  the  definition  of  "compensating  ratio".   The 

constants,  b.  and  c.  are  determined  from  the  coefficients  of 

3  1 

the  characteristic  equation.   Since  for  this  case,  p.  is  real, 
the  compensating  ratio  is: 

-TE-  =  -Constant  (5-13) 

Example  5-1 

Figure  5-4  is  a  block  diagram  of  a  type  one,  unity  feed- 
back control  system.   The  system's  characteristic  equation  is 
third  order  and  contains  the  variable  parameters  a  and  3. 
The  characteristic  equation  of  the  system  shown  in  figure  5-4 


is  : 


s3  +  (P1+8)s2   +  (Px3+K)s  +  aK  =  0         (5-14) 
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COMPENSATION 


PLANT 


Figure  5-4:   Block  Diagram  Example  5-1  System 

Choosing  the  components  and  gain  factor  such  that  the  system 
is  stable  let: 

P,  =  12;  a    =  3;  3  =  1;  K  =  210. 

The  characteristic  equation  is  now: 

s3  +  13s2  +  222s  +  630  =  0 
which  factors  to: 

(s+3.318)  (S+4.84  +  J12.9)  (s+4 . 84- j 12 . 9)        (5-15) 
Without  any  type  of  compensation  or  self-adaption,  the 
behavior  of  the  system  is  dependent  upon  changes  in  the  param- 
eter 3.   Table  5-1  lists  various  values  of  3  and  the  corres- 
ponding factors  of  the  systems  characteristic  equation.   For 
convenience  these  values  are  plotted  in  figure  5-5. 

Table  5-1:   System  Roots  With  a  Variable  Parameter 

3  Values  Factors  of  Characteristic  Equation 

.01  (s+3.49  )  (s+4.26+jl2.73)  (s+4 . 26- jl2 . 73) 

.1  (s+3.47  )  (s+4.31+jl2.75)  (s+4 . 31- J12 . 75) 

1.0  (s+3.318)  (S+4.84+J12.9  )  (s+4 . 84- j 12 . 9  ) 

2.0  (s+3.1   )  (s+5.4  +J13.0  )  (s+5.4  -J13.0  ) 

10.0  (s+2.2   )  (s+9.9  +J13.7  )  (s+9.9  -J13.7  ) 
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Figure  5-5:   System  Poles  With  a  Variable  Parameter 


Assuming  the  desired  operating  point  for  the  system  is 
that  position  specified  by  the  value  3  =  1^   equation  5-15 
gives  the  factors  of  the  characteristic  equation.   Note  that 
the  dominant  root  of  the  system  is  p.  =  -3.318.   The  compensa- 
tion network  has  been  added  to  provide  adaption  so  that  any 
variations  in  3  can  be  cancelled  by  the  adjustment  of  the 
parameter  a  insofar  as  the  dominant  root  performance  is  con- 
cerned. 

Without  specifying  any  change  in  a  or  3 ,  evaluate  the 
constant  term  of  equation  5-13,  i.e., 


2 .         n 

cQ+clPi+c2Pi  +...+cnPi 

b0  +  b1pi+b2Pi24-...+bnPin 


(5-17) 
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aQ  =  Ka 


bQ  =  K  =  210;  cQ  =  dQ 


a1   =  +  P,0  +  K 


b,  =  0;  c,  =  12;  d,  =  K  =  210 


a2  =    +  3  +  P-l 


a3  -  1 


b2  =  0;  c2  =  1;  d2  =  P   =  12 


b3  =  c3  =  0;  d3  =  1 


With  the  value  p.  =  -3.318,  equation  5-12  becomes: 

da  =   0+12(-3.318)+l(-3.318)2+0 
d~3  "         210+0 

reducing 

da     28.7936 
d3       210 

hence : 

da  =  +.13711d6  (5-18) 

From  the  results  of  equation  5-18,  for  any  change  in  the 
plant's  parameter  3/  in  order  to  keep  the  root  p.  =  -3.318 
constant,  a  must  be  changed  (.137)d$. 
Situation  (1) : 

Assume  3changes  from  a  value  of  1  to  a  new  value,  3=2. 
As  the  change  in  3  is  plus  one,  a  must  be  changed  (.137)  (1) 
to  a  value  a1  =  +3.137.   These  new  values  of  a  and  3  change 
the  characteristic  equation  to: 

s3  +  14s2  +  234s  +  658.79  =  0 
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By  a  root  solving  method,  the  above  equation  factors  to: 

(s+3.3179) (s+5.34+jl3.04) (s+5 . 34- J13. 04) 
The  dominant  root  has  been  held  constant. 
Situation  (2) : 

Assume  $  changes  from  a  value  of  1  to  a  value  6  =  10. 
The  variable  parameter  a  must  then  be  changed  (.137)  (9)  = 
1.233  to  a  new  value  a'  =4.233.   The  new  characteristic  equa- 
tion becomes : 

s3  +  22s2  +  330s  +  889.14  =  0. 
This  equation  factors  to: 

(s+3.3174) (S+9.34+J13.45) (s+9 . 34- jl3. 45) 
The  dominant  root  has  been  held  constant. 
Situation  (3) : 

Assume  $  changes  from  a  value  of  1  to  $  =  .9.   This  10 
percent  decrease  forces  a  to  change  -.0137  to  a  new  value 
a'  =  2.98  73.   The  characteristic  equation  becomes: 

s3  +  12.9s2  +  220.8s  +  627.123  =  0 
which  factors  to: 

(s+3.318) (s+4.79+jl2.88) (s+4 . 79- jl2 . 88) . 
The  dominant  root  has  been  held  constant. 
Example  5-2 

Given  a  type  '0',  unity  feedback  system  as  shown  in  fig- 
ure 5-6.   The  system  characteristic  equation  is  third  order 
and  is  given  by: 

S3  +  (P1+P2+3)s2  +  [P1P2+KK1+(P1+P2)  g]s  +  KKja  +  P-^S  =  ° 

(5-19) 


73 


R(s) 


^             <s 

K1 (s+a) 

t= < — 

(s+P;L) 
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K 
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^ 

(s+P2)  (s+3) 

Figure  5-6:   Block  Diagram  Example  5-2  System. 

Let  the  values  of  the  various  components  be  assigned  such  that 
the  system  is  stable.   These  values  are: 

P1  =  2;  P2  =  3;  KK-  =  4;  a  =  5;  3  =  6. 

The  characteristic  equation  with  these  values  substituted  be- 
comes : 

s3  +  lis2  +  40s  +  56  =  0. 
The  factors  of  the  above  equation  are: 

(s+5.716)  (s+2.642+jl.678)  (s+2 . 642-jl. 678) . 
The  results  of  the  Case  (1)  compensating  ratio  can  be 
used  such  that  for  any  variation  in  the  parameter  3,  one  will 
be  able  to  vary  the  parameter  a  and  keep  the  real  root,  d.  = 
-5.716  constant.   The  procedure  to  be  followed  can  be  outlined 
into  three  steps. 

Step  (1) :   Determine  the  value  of  the  constant  terms  in  the 
coefficients  of  the  characteristic  equation: 
aQ  =  KK^^a  +  P!P2e 

bQ  =  KK1  =  4;  cQ  =  PlP2  =  6;  dQ  =  0 
al  =  P1P2  +  KK1  +  (P1+P2)3 


b,  =  0;  c,  =  P,  +  P0  =  5;  d,  =  P,P0  +  KK.  =  10 


V  2 


(5-20) 
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a2  =  Px  +  P2  +  3 


a3  =  1 


b2  =  0;  c2  -  1;  d2  =  P;L  +  P2  =  5 


b3  =  c3  *•  0;  d3  =  1  (5-20) 


Step  (2) :   Using  the  values  found  in  equation  5-20  and  the 
real  root  p.  =  -5.716,  evaluate  the  compensating  ratio  equa- 
tion 5-12. 

da     6+5(-5.716)+l(-5.716)2 

d3  ~  4" 

||  =  -9'°^26  or  da  =  -2.273dB  (5-21) 

Step  (3):   Identify  the  change  in  the  parameter  3  and  adjust 
the  value  of  a  according  to  equation  5-21.   The  real  root 
will  remain  constant. 

Assume  3  changes  from  a  value  of  +6  to  a  new  value  of 
5.3.   d3  =  -.7  and  the  change  in  a  must  be: 

da  =  -2.273(-.7)  =  1.59 
Substitute  the  new  value  of  a  =  5  +  1.59  into  the  character- 
istic equation  along  with  the  new  value  of  3.   With  a  =  6.59 
and  3  =  5.3,  the  new  characteristic  equation  becomes: 

s3  +  10.3s2  +  36.5  +  58.2  =  0. 
This  equation  factors  to: 

(s+5.693)  (s+2.31+j2.214)  (s+2 . 31-J2 . 214)  . 
5.4:  Case  (2) :   The  Dominant  Roots  are  a  Complex  Conjugate 
Pair  and  are  to  be  Held  Constant 

When  the  dominant  roots  of  the  characteristic  equation 
are  a  complex  conjugate  pair,  an  increment  in  one  parameter, 
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Figure  5-7.   Change  in  Root  Location  With  a 
Change  in  a  System  Parameter 


such  as  ft,  changes  the  root  locations  as  shown  in  figure  5-7. 

Although  any  new  root  location  may  occur,  it  is  impor- 
tant to  note  that  the  two  roots  will  always  be  complex  con- 
jugates.  The  behavior  of  complex  conjugate  roots  is  far  more 
complex  than  the  real  root  case  presented  earlier.   Instead 
of  a  simple  plus  or  minus  variation  as  in  the  real  root  loca- 
tion, an  uncontrollable  parameter  variation  in  the  system's 
plant  will  cause  the  complex  dominant  roots,  p.  and  p.*,  to 
change  both  in  their  real  and  imaginary  cormoonents.   In  this 
section,  it  is  proven  that  with  two  variable  parameters  in 
a  system,  the  dominant  root  locations  can  be  returned  to 
their  original  positions  by  corrective  compensation  in  only 
very  trivial  cases. 

For  convenience  equation  5-11  is  repeated  below.   This 

equation  shows  that  the  differential  of  any  root,  p.,  is  a 

function  of  the  change  in  any  two  variable  parameters  a   and 

B. 
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b0+blP.-K..+bnp.n  C0+C1p.+...+CnPin 

dPi  =  ^Tp-1 da  +   h^TT d6 

(5-11) 
where:   b.  and  c.,  j  =  0,1,. ..,11  are  constants  determined  by 
the  system, 

:   a  and  3  are  real  parameters,  initially  fixed  in 
value  but  subject  to  variations. 

If  p.  is  complex  and  is  desired  to  be  held  constant, 
then  the  differential  of  the  root  p.  must  be  set  equal  to 
zero ,  i.e., 

dPi  =  0  . 

Equation  4-11  then  reduces  to: 

da    co+clPi+-"+cnpin 

S|  =  — — S-i-  (5-22) 

at5   bn+blP.  +  .  .  ,+b  p.n 
0   l^i      n*i 

As  the  parameters  a  and  3  are  specified  to  be  real  by 
definition,  then  the  right-hand  side  of  equation  5-22  must 
be  real  to  have  any  meaning.   There  must  be  restrictions  on 
the  system  parameters  in  order  to  make  the  compensating  ratio, 
equation  5-22,  equal  to  a  real  constant. 

To  find  these  restrictions,  assume  the  numerator  of  the 

compensating  ratio  with  a  complex  root,  p. ,  may  be  combined 

into  the  general  format:   m  +  jn.   Likewise  assume  that  the 

denominator  is  expressable  as:   x  +  jy.   Equation  5-21  mav 

now  be  written: 

da  m   _m+4n  (5_23) 

d3     x+jy 
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Separating  the  above  equation  into  its  real  and  imaginary 

components : 

da     (m+jn)  .  (x-jy)  (5_24) 

dT     (x+jy)    (x-jy)  ls  *qi 


or 


da      mx+ny    . nx-my 
d3  =     x2--     -'-2... 


(5-25) 


x  +y      x  +y 

Equation  5-25  will  be  real  if  and  only  if  nx  =  my.   The 

compensating  ratio  is : 

da    m    n  ,c    ~>r\ 

M  =   "x  =  ~y  (5"26) 

The  numerator  and  the  denominator  of  the  compensating  ratio 
must  be  constant  factors  of  each  other.   As  will  be  shown, 
when  p.  is  complex  the  above  conditions  imply  pole-zero  can- 
cellation.  This  result  is  expressed  in  the  following  theorem, 
THEOREM  5-1 

—  "f-Vi 

For  the  n    order  feedback  control  system  shown  in  fig- 
ure 5-8  with  variable  system  parameters,  a  and  3,  there  will 
exist  a  closed-loop  root,  p.,  whose  location  remains  invar- 
ient  for  all  increments  in  a  or  3  such  that  the  compensating 
ratio,  da/d3,  is  equal  to  a  real  constant  provided,  if  and 
only  if: 

(a)  the  root  of  interest,  p. ,  is  real  or 

(b)  the  variable  parameter  a  equals  the  variable 
parameter  3  which  implies  a  zero-pole  cancel- 
lation, a  trivial  case. 
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Figure  5-8.   n    Order  Feedback  Control  System 


Proof: 


For  convenience,  the  compensating  ratio  for  any  n 


-th 


order  feedback  control  system  is  repeated: 

n 


da 
d(3 


cn+c, p. +. . ,+c  p 
0   l-i n- i 

n 


(5-22) 


bA+b,p. +. . .+b  p 
0   1*1      n  i 

If  the  above  ratio  is  to  be  equal  to  a  real  constant  'K' , 

then  equation  5-22  reduces  to: 


K(b0+blP;L+.  .  -+bnPin)  =  -(c0+clPi+.  ..+cnDin) 


equating  to  zero: 


n 


(bnK+cn)  +  (b.K+cJp.  +  ...  +  (b  K+c  )p.   =  0 
0    0        1 


(5-27a) 


(5-27b) 


l'1"!  '  n  ~n'  •'  l 

A  sufficient  condition  for  part  (a)  of  Theorem  5-1  is 
that  equation  5-27  is  always  satisfied  when  the  root  of  in- 
terest, p. ,  is  real.   This  can  be  shown  to  be  true  by  insoec- 
tion  of  the  compensating  ratio,  equation  5-22.   The  param- 
eters a  and  $  are  by  definition,  real  parameters.   The  con- 
stant terms,  b.  and  c.,  are  also  defined  to  be  real  numbers. 
If  p.  is  real,  then  equation  5-22  is  a  ratio  of  two  real 
numbers,  hence  da/d$  is  equal  to  a  real  constant.   If  this 
real  constant  is  equal  to  'K1,  then  equation  5-27  must  alwavs 
be  satisfied  for  any  real  root,  p.. 
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A  sufficient  condition  for  part  (b)  of  Theorem  4-1  fol- 
lows since  if  a  =  B,  then 

s+a 


s+B 


=  K  =  1 


(5-28) 


The  closed  loop  system  in  figure  5-8  is  independent  of  a  and 

b. 

The  necessary  condition  for  the  proof  of  Theorem  5-1  can 
be  shown  with  the  aid  of  figure  5-9  where  the  real  factor,  k, 
has  been  added  for  generality. 


R(s)  t 


"N     „ 

G1(s) 

\  G2(s)  -  m  k 

C(x) 

J     * 

? 

r 

Figure  5-9.   Modified  Closed  Loop  Feedback  Control  System 


The  closed-loop  control  system  of  figure  5-8  has  been 
transformed  into  the  modified  control  system  of  figure  5-9 
by  the  following  modifications: 

G,(s)  =  G(s)/H(s),  the  fixed  portion  of  the  closed 
loop  system. 

G2(s)  =  k (s+a) / (s+B) ,  the  adjustable  portion  of  the 
closed  loop  system. 
The  characteristic  equation  is: 


or 


or 


k  |(s)  (s+o)  +   m 
H(s)  (s+B) 


H(s)  (s+B)  +  kG(s)  (s+a)  =  0 


s[kG(s)+H(s) ]  +  akG(s)  +  BH(s)  =  0 
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(5-29a) 


(5-29b) 


(5-29c) 


Assume  that  G(s)  and  H(s)  may  be  expressed  as: 

G(s)  =  gQ  +  g,s  +  ...  +  g  _,s  ~  (5-30a) 


and 


H(s)  =  hn  +  h,s  +  ...  +  h   ,sn  1  (5-30b) 

u    1  n-i 


Substituting  equations  5-30a  and  5-30b  into  equation  5-29c, 
the  characteristic  equation  becomes: 

(kg0a+hQ3)  +  (kg^+h^S+kgQ+hQ)  s  +  ... 

+  ^B-A-l^B-2+V2,|aJi 

+  (kgM  ,+h   ,)sn  (5-31) 

n-l  n-l 

— +■  v» 
Equation  5-31  is  the  n    order  characteristic  equation  of 

the  feedback  control  system  of  figure  5-8,  which  is  written 

as  : 

an   +   a,s  +  ...  +  a   , sn   +  a  sn  =  0       (5-32) 
o    1  n-l        n 

From  the  introduction,  the  coefficients  of  the  characteristic 

equation  are  expressable  as: 

a.  =  b.a  +  c .8  +  d.  (5-4) 

J     D      3  D 

j  =  0,1,. . . ,n 

Equating  the  coefficients  of  equations  5-31  and  5-32  yields: 

aQ   =  kgQa  +  hQe  =  bQa  +  CqB  +  dQ   (5-32a) 

al    =  k^ia  +  hi^  +  ^k^n+h0^  =  bla  +  cl^  +  dl   (5-32b) 

•    •    • 

a   ,  =  kg   ,  a  +  h   ,  6  +  (kg   0+h   0)=b   ,a  +  c   ..  6  +  d   , 
n-l    ^n-1     n-l      ^n-2   n-2     n-l     n-l     n-l 

(5-32c) 

a^   =  (kgn   +h   r)  *  b  a  +  c  B  +  dn   (5-32d) 

n        n-l   n-l  n     n     n 


81 


Equating  the  constant  terms : 


=  kg. 


=  kg. 


=  h" 


=  0 


di   =  kqo  +  ho 


(5-34a) 
(5-34b) 


b   t  =  kg 


n=l 


n 


n-1 


=  0 


c   n  =  h   n 
n-1     n-1 


n 


=  0 


dn-l  =  k*n-2  +  hn-2(5"34c) 


I    =  kg   ,  +  h   , (5-34d) 
n      3n-l    n-1 


With  the  above  values  of  b.  and  c..  equation  5-27b  becomes: 

D       3  r 

(kgQK+h^  +  (kg1K+h1)pi  +  ...  +  (kgn_1K+hn_1) pin"1  =  0     (5-35) 


*® 


G(s) 
HTiT 


K 


Figure  5-10.   (n-1)    Order  Feedback  Control  System 

Figure  5-10  shows  an  (n-1)    order  closed  loop  feedback 

control  system.   The  characteristic  equation  of  this  system 

is : 

K  ir^T  +1=0  (5-36a) 

H  (s) 

or 

KG(s)  +  H(s)  =  0  (5-36b) 

If  G(s)  and  H(s)  are  represented  by  equations  5-30,  then 
equation  5-36b  can  be  written  as: 

(g0K+hQ)  +  (g1K+h1)s  +  ...  +  (g^^"^..-^  s11"1  =  0  (5-37) 
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Comparing  equation  5-35  with  equation  5-37,  it  follows  that 
the  root,  p. ,  must  also  be  a  root  of  the  system  shown  in  fig- 
ure 5-10,  that  is, 

kK  !j|j.  +  1  =  0  (5-38) 

Since  p.  is  a  root  of  the  original  system  shown  in  figure 
5-8,  then  it  also  follows  that: 

T?WlTfT+1  =  0  (5"39) 


1 

Equating  equations  5-38  and  5-39: 


(p.+a) 

TiTW  =  kK  (5-40) 


Solving  for  the  value  of  the  root,  p. : 


(p.+a)  =  kK(p.+3)  (5-41) 


or 

0-a 


i   1-kK 


(5-42) 


From  equation  5-42,  the  root  of  interest,  p.,  is  always 
real  if  kK  is  not  equal  to  one.   If  kK  is  equal  to  one,  then 
p.  is  indeterminate,  however  from  equation  5-28,  kK  =  1  im- 
plies a  =  3  or  zero-pole  cancellation,  a  trivial  case  of 
compensation . 

5.5:  Case  (3):   The  Dominant  Roots  are  Complex  Conjugates  and 

the  Damping  Coefficient  is  to  be  Held  Constant 

Since  the  damping  coefficient,  cu)  ,  is  to  be  held  con- 
stant for  variations  in  the  plant's  parameter  3/  let  the  real 
part  of  either  dominant  complex  conjugate  root  be  held  con- 
stant, i.e., 


<£ 


[p. ]  =  CONSTANT  (5-4  3) 

e  ltiJ 
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Q. 


With  the  If/    [p.]  equal  to  a  constant  value  the  real  part  of 


e  ■  1 


the  differential  of  p.  must  equal  zero,  or 

^e[dp.]  =  0  (5-44) 

With  the  above  condition,  the  real  part  of  equation  5-11  is 
written: 


ftw-«-g.[v*q;ij*' 

^1 VW-ii 


da 


dB 


(5-45) 


hence : 


«■ 


bQ+blPi+---+bnPi 
-F'Cp.) 


n 


da 


-R, 


c~+c, p.+. . ,+c  p. 
0   1"  l n-  i 


n 


d$ 


(5-46) 


Evaluating  the  terms  in  the  parenthesis  of  equation  5-46,  let 

b„+b,p.+. . ,+b  p.n 
m  +  jn  =   °   i£k-  v   n  x  (5-47) 


and 


u  +  jv 


-F'(Pi> 


c.+c, p.+. . .+c  p. 
0   l-i      n-  l 


n 


-F'(p-) 


(5-48) 


Using  the  above  notation,  equation  5-16  is  rewritten  as: 

/Oe[mda+jnda]  =  -(/JQ  [ud3+jvd6]  (5-49a) 
or 

mda  =  ud$  (5-49b) 

The  compensating  ratio  for  a  constant  damping  coefficient  is 
therefore  equal  to: 

at  "  -i  '5-50» 
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Example  5-3 

Using  the  system  of  example  5-2,  the  characteristic  equa- 
tion is: 

s3  +  (P1+P2+3)s2  +  [P1P2+KK1+(P1+P2)  3]  s  +  Kl^a  +  P-^23  =  0 

with  the  parameter  values: 

P1  =  2;  P2  =  3;  KK1  =  4;  a  =  5;  6  =  6. 

The  numerical  characteristic  equation  is 

s3  +  lis2  +  40s  +  56  =  0 
which  factors  to: 

(s+5.716)  (s+2.642  +  jl.678)  (s+2 . 642-jl. 678)  . 
The  parameter  3  will  vary  and  again  must  be  compensated 
for  by  a  change  in  the  parameter  a,  such  that  the  damping 
coefficient,  £uj  ,  is  held  constant.   The  damping  coefficient 
will  be  held  constant  if  equation  5-50  is  satisfied,  there- 
fore the  following  procedure  is  recommended. 

Step  (1) :   Determine  the  value  of  the  constant  terms  in  the 
coefficients  of  the  characteristic  equation. 

This  step  has  been  performed  in  example  5-2  and  equa- 
tion 5-20  is  repeated  for  convenience. 

bQ  =  4      cQ  =  6     dQ  =   0 

b,  =  0     a,  ■  5     d1  =  10 

b2  =  0      c2  =  1     d2  =   5 

b3  =  0      c3  =  0     d3  =   1 

Step  (2) :   Form  the  real  part  of  equations  5-47  and  5-48. 


Denominator: 

3       2 
F(s)    =  s   +  lis   +  40s  +  56 

F*  (s)   =  3s2  +  22s  +  40 

p.      =  -2.64  -  jlO.35 

F'(p.)  =  3(-2.64-jl.68)2  +  22  (-2 . 64- jl . 68)  +  40 
F' (p. )  =  -5.64  -  jlO.35 

Equation  5-47:  2 

b0+blPi+blP. 

m  +  ^n  =  ^Tp-1 

4 
~  5.64+J10.35 

m     ■             4(5.64)  les 

m  =  s— 7T  m  =  .163 

(5.64r+(10.35r 

Equation  5-4  8:  2 

C-.+c.p.+c-p . 
0   l^i   2^i 

U+  ^V=  -F'(p.) 


-3.1+1.43  _   3.13  /1_^ 
5.64+J10.35  "  11.8  /61.T 

u  =  -.094 

The  compensating  ratio  is  the  negative  ratio  of  u  over  m  or 

da  =  +.573d$ 
Step  (3) :   Identify  the  change  in  the  parameter  3  and  adjust 

Assume  3  changes  from  a  value  of  6  to  a  value  of  5.3. 
The  change  in  3  is: 

d3  =  -.7 

da  =  (-.7)  (.573) 

da  =  -.4011 
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The  parameter  a   must  change  to  5  -  .4011  =  4.59  8?,  The  new 
characteristic  equation  will  be: 

s3  +  10.3s2  +  36.5s  +  50.2  =  0. 
This  equation  factors  to: 

(s+5.01) (s+2.64+jl.74) (s+2 . 64- j 1 . 74) . 
By  inspection,  the  real  part  of  the  complex  conjugate  pair 
has  been  kept  constant. 

5.6:  Case  (4)  :   The  Dominant  Roots  are  a  Complex  Con-jugate 

Pair  and  the  Damped  Resonant  Frequency  is 
to  be  Held  Constant 
In  a  similar  manner  as  the  damping  coefficient,  £w  ,  was 
held  constant  in  case  (3),  the  damped  resonant  or  natural 
frequence,  oo , ,  will  be  held  constant  by  keeping  the  imaginary 
part  of  either  dominant  complex  root  constant,  i.e., 

<-j]  [p. ]  =  CONSTANT 

With  the  Imaginary  part  of  p.  held  constant,  the  imaginary 
part  of  the  differential  of  p.  must  equal  zero,  or 


4 


fdp.]  =  0  (5-51) 

m  c  l 


Setting  the  imaginary  part  of  the  differential  of  p.  equal 

to  zero,  equation  5-11  becomes: 

b-.+b,  p.  +  . . .+b  p . 
0   l*i      n^i 


4 


in 


■F'(Pi) 


n 

da 


■4 


c0+c1P.+...+cnPin 

=FT?7! de 


m 

i_  ""  i 

(5-52) 
Inspecting  the  terms  in  the  parenthesis  as  done  in  case  (3) , 
equations  5-47  and  5-48  are  evaluated.   The  comDensatina 
ratio  for  a  constant  damped  resonant  frequency  is  found  to  be 

da         V  ,r   r,v 

dT=  'n  (5"53) 
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where 


n 


and 


A    c0+clPi+...+cnp. 
V  =  ^m   -F'(Pj) 


A   fbQ+blPi+...+bnPin 
n  _^m  L -F'  (p.) _ 


(5-54) 


(5-55) 


Example  5-4 

Using  the  system  of  example  5-2  and  5-3,  the  characteris- 
tic equation  is : 

s3  +  (P1+P2+6)s2  +  [P;,P2+KK1+(P1+P2)  Bis  +  KKxa  +  P-^S  =  0 

Substituting  in  the  values  of  the  parameters  as  before,  the 
characteristic  equation  becomes : 

s3  +  lis2  +  40s  +  56  =  0 
which  factors  to: 

(s+5.716)  (s+2.64+jl.68)  (s+2. 64-jl. 68) . 
In  this  example  the  damped  resonant  frequency  will  be 
held  constant  by  the  following  procedure  as  indicated  in  the 
above  case. 

Step  (1) :   Determine  the  value  of  the  constant  terms  in  the 
coefficients  of  the  characteristic  equation. 
Repeating  equation  5-20  for  convenience: 


bQ  =  4 


b1  =  0 


cQ  =  6 


c,  =  5 


dn   =   0 
0 


di  =  10 


b2  =  0 


c2  ■  1 


b3  =  0 


c3  =  0 


d3=   1 


Step  (2) :   Form  the  imaginary  parts  of  equation  5-47. 
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Denominator 


F(s)   =  s3  +  lis2  +  40s  +  56 

F'  (s)  =  3s2  +  22s  +  40 

p.     =  -2 . 64-jlO . 35 :  the  dominant  root  of 


interest.   As  before  in  example  5-2: 


Equation  5-47: 


F' (p.)  =  -5.64-jl0.35. 


m  +  jn  =  5.64+J10.35 


4 
11.8/61.5" 


n  =  -.298 


Equation  5-48: 


-3.1+J.43 
U  +  JV  ~    5.*4+;jl0.35 

=  .266/110o46' 

v  =  .246 

Step  (3):   Determine  the  compensating  ratio  according  to 

equation  5-53. 

da      .246 
d3    "-.298 

da  =  +.825d3 

Step  (4):   Identify  the  change  in  the  parameter  3  and  adjust 

If  3  changes  from  a  value  of  6  to  a  new  value  3  =  5.3, 

the  change  in  3  is: 

d3  =  -.7 

da  =  (-.7)  (.825) 

da  =  -.577 
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The  parameter  a  must  change  to 

5  -  .577  =  4.423  . 
The  new  characteristic  equation  will  be: 

s3  +  10.3s2  +  36.5s  +  49.4  =  0. 
This  equation  factors  to: 

(s+4.94)  (s+2.68+jl.683)  (s+2 .68-jl. 683)  . 
By  inspection  the  imaginary  part  of  the  dominant  root 
has  been  held  constant  by  proper  use  of  the  compensating  ra- 
tio, equation  5-53. 
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CHAPTER  VI 
Conclusions 

In  this  thesis,  a  purely  mathematical  property  of  all 
polynomials  of  degree  n,  (called  the  associated  polynomials 
of  the  defining  polynomial)  has  been  applied  to  the  fields 
of  circuits  and  systems.   The  properties  of  associated  polyno- 
mials have  been  investigated  and  examples  given  that  demon- 
strate that  these  polynomials  provide  useful  results.   It 
has  been  shown  that  associated  polynomials  may  be  used  in  a 
new  method  for  accomplishing  a  partial  fraction  expansion  in- 
cluding the  troublesome  case  of  repeated  roots  of  high  order 
of  multiplicity.   Associated  polynomial  theory  has  been  used 
to  study  the  relationship  between  variations  of  polynomial 
coefficients  and  polynomial  roots.   The  theorems  developed 
in  Chapter  Three  that  show  that  the  variation  of  the  coeffi- 
cients of  a  polynomial  with  respect  to  the  variation  of  its 
roots  or  vice-versa  is  determined  by  associated  polynomials 
have  important  engineering  applications.   They  are  immedi- 
ately used  in  a  root  solving  procedure  and  to  express  sensi- 
tivity coefficients  in  a  new  analytical  form. 

The  most  significant  result   of  this  thesis  is  the 
investigation  of  the  theory  of  compensating  parameters  as 
applied  to  self-adaptive  control  systems.   A  mathematical 
relationship  between  two  variable  parameters^ such  that  the 
system's  dominant  root  performance,  has  been  derived.   In 
addition,  a  theorem  of  fundamental  importance  in  adaptive 
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control  theory  has  been  presented  that  Effects  the  basic 
theory  of  complex- con jug ate  dominant  root  adaptive  compensa- 
tion. 

As  a  suggestion  for  further  research,  it  is  recommended 
that  an  actual  computer  program  be  written  that  performs  the 
partial  fraction  expansion  of  any  ratio  of  polynomials.   In 
addition,  it  is  believed  that  the  results  of  Chapter  Three 
have  logical  extensions  to  a  derivation  of  an  analytical  ex- 
pression for  a  root-locus  plot.   This  analytical  expression 
could  then  be  programmed  onto  a  digital  computer.   Lastly, 
the  results  of  Chapter  Five  can  be  extended  to  variable 
parameters  that  are  not  linearly  related  to  the  coefficients 
of  the  characteristic  equation,  i.e.,  the  variable  parameters 
would  appear  in  the  form: 

a.  =  b.a  +  c.3  +  h.a$  +  d.  (6-1) 

113  3  3 


or: 


2  2 

a.  =  b.~a   +  b..a  +  h.a$  +  c.,$  +  c.~3   +  d.    (6-2) 
3  :2       jl  j  jl  -\2  j 
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APPENDIX  A 
Horner's  Method 
A .  1 :  Introduction 

A  complete  historical  survey  of  associated  polynomials 
is  given  in  reference  1.   In  this  reference  it  is  noted  that 
what  has  been  labeled  associated  polynomials  are  sometimes 
referred  to  as  Horner  polynomials.   This  appellation  comes 

about  by  three  processes  involving  associated  polynomials 

(4) 
that  are  attributed  to  Horner    .   These  processes  are: 


(1) 
(2) 
(3) 


The  nested  evaluation  of  polynomials; 
Synthetic  division; 

A  method  of  approximating  the  roots  of  polyno- 
mial equations. 


Each  of  the  three  processes  above  are  commonly  called  Horner's 
method  or  rule     .   As  is  shown,  the  method  of  approximating 
the  roots  of  a  polynomial  directly  involves  synthetic  divi- 
sion. 
A. 2:  The  Nested  Evaluation  of  Polynomials 

Associated  polynomials  may  be  used  in  the  evaluation  of 
polynomials  by  a  process  called  "nesting".   In  order  to  eval- 
uate the  polynomial, 

F(s)  =  a  s   +  a   ,s     +  ...  +  an :   a   =  1     (1-1) 
n      n-1  0     n 

for  any  value  of  s,  there  is  required  2n  +  1  multiplications 
and  n  additions.   Using  associated  polynomials  there  is  re- 
quired only  n  multiplications  and  n  additions.   This  process 
of  evaluation  is  performed  by  the  following  method: 

F(s=u)  =  a  un  +  a   ,un~   +  . . .  +  an :   a   =1    (A-I) 
n      n-l  u    n 
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but 


where : 


F(u)  =  AQ(u)  (1-12) 


AQ(u)    =  uA-^u)  +  aQ  (l-9n) 

Ax(u)    =  uA2(u)  +  a,  (l-9n-l) 


A^  . (u)  =  uA  (u)  +  an    ,  (l-9b) 

n-1        n       n-1 


and  since: 


An+1(u)  =  0  (1-7) 

A  (u)    =  uA.  ,  (u)  +  an   =  aM  (l-9a) 

n  n+i       n    n 

Hence  F(u)  is  evaluated  using  only  n  multiplications  and  n 

additions . 

A.  3:  Synthetic  Division 

The  equation: 

_,  N    An(u)     n         .  , 

^1=  -2 +  I      A^uJs^1  (1-15) 

s-u      s-u     •  -i   1 

3  =  1   J 

states  that  the  quotient  and  remainder  after  division  by  a 

polynomial  of  the  form  s-u  may  be  calculated  recursively. 

The  nature  of  the  recursion  is  such  that  one  may  work  with 

detached  coefficients.   By  detached  coefficients  is  meant  an 

abbreviation  of  the  ordinary  multiplication  and  division 

processes  used  in  ordinary  algebra.   The  coefficients  alone 

with  their  signs  are  used  in  the  process.   The  powers  of  the 

variable  occurring  in  the  various  terms  are  understood  only 

from  the  order  in  which  the  coefficients  are  written.   Missing 
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powers  are  assumed  to  be  present  with  zero  coefficients.   The 
process  of  synthetic  division  is  best  understood  by  the  follow- 
ing example  and  summary. 
Example  A-l 
Given  the  polynomial: 

3      2 
F(s)  =  As   +  BS   +  Cs  +  D  (A-2) 

And  the  divisor  polynomial: 

s  -  a  (A-3) 

The  usual  process  of  division  is  as  follows: 

3      2 

As   +  Bs   +  Cs  +  D/   s-a 

3       2  2 

As   -  Aas  As  +Es+F 

2 

Es   +  Cs 

2 

Es   -  Eas 


Fs  +  D 
Fs  -  Fa 


Or: 


D  +  Fa  =  R 


AS3+BS2+CS+D  _  (As2+Es+F)  +  _J_  (A.4) 


s-a 
Comparing  with  equation  1-15: 


and 


AQ(a)  =  R  (A-5 


A.UJs3    =  (As  +Es+F) :   j  =  1,2,...  ,n       (A-6) 


Using  synthetic  division,  the  process  of  division  is  reduced 

by  only  using  the  coefficients  according  to  the  following 

scheme. 

A   B    C    D   /+a 

Aa   Ea   Fa 
A   E    F    R 
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Notice  that  the  first  three  numbers  of  the  last  line, 

2 

A,  E,  and  F,  are  the  coefficients  of  s  ,  s,  and  1,  in  the  Quo- 
tient, and  R  is  the  remainder.   The  process  may  be  expanded 
to  any  n    order  polynomial. 

As  a  summary,  synthetic  division  may  be  performed  accord- 
ing to  the  following  rule. 

(1) :   To  divide  F(s)  by  a  polynomial  of  the  form,  s  -  a, 
arrange  F(s)  into  descending  powers  of  s. 

(2):   Write  the  coefficients  of  F(s)  on  a  horizontal 
line,  in  the  order  which  corresponds  to  the  arrangement  spec- 
ified in  (1) .   If  any  power  of  s  is  missing  in  F(s) ,  supply 
that  power  with  a  zero  coefficient. 

(3) :   Multiply  the  first  coefficient  A  by  a,  write  the 
product  below  the  second  coefficient  B  and  add.   Multiply 
this  sum,  E,  by  a,  write  the  product  below  the  third  coeffi- 
cient, D,  and  add. 

(4) :   Proceed  in  this  way  until  all  of  the  places  in  the 
third  row  except  the  first  are  filled  up,  and  write  down  the 
first  coefficient,  A,  of  F(s)  in  the  first  place  of  the  third 
row.   The  last  number  of  the  third  row  will  be  the  remainder, 
and  the  other  numbers  in  the  third  row  will  be  the  coeffi- 
cients of  the  quotient  obtained  when  F(s)  is  divided  by  s  -  a. 
A. 4:  Approximation  of  the  Roots  of  a  Polynomial  Equation 

A  method  for  approximating  the  real  roots  of  an  algebraic 

equation  first  appeared  in  a  paper  read  by  W.  G.  Horner  be- 

(4) 
fore  the  Royal  Society  in  1819    .   This  process  involves 

associated  polynomials  in  that  repeated  use  of  synthetic 
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division  is  used  to  expedite  the  work  of  substituting  values 

(21) 
for  the  variable.   The  essential  steps     of  Horner's  method 

of  approximation  are  as  follows: 

(1):   Using  curve  plotting,  the  theorem  of  continuous 
fractions*/  etc.,  isolate  a  positive  root  between  two  succes- 
sive integers.   If  the  equation  has  only  negative  roots, 
transform  it  to  one  whose  roots  are  the  negatives  of  those 
of  the  given  equation. 

(2) :   Transform  the  equation  into  an  equation  whose  roots 
are  decreased  by  the  lesser  of  the  integers  between  which  the 
root  lies,  by  the  substitution  x'  =  x-a.   The  root  of  the  new 
equation  will  lie  between  zero  and  unity. 

(3):   Isolate  the  root  of  the  new  equation  between 
successive  tenths. 

(4) :   Transform  the  last  equation  into  an  equation  whose 
roots  are  decreased  by  the  smaller  of  these  tenths  and  iso- 
late the  root  of  this  equation  between  hundredths.   Continue 
this  process  to  one  decimal  place  more  than  the  place  to 
which  the  answer  is  to  be  correct. 

The  root  sought  is  then  the  total  amount  by  which  the 
roots  of  the  original  equation  were  reduced,  namely,  the 
lesser  integer,  plus  the  lesser  tenth,  plus  the  lesser  hun- 
dredth, etc.,  the  last  decimal  being  rounded  off  to  make  the 
result  accurate  to  the  desired  decimal  place. 


*If  P(x)  is  a  polynomial  with  real  coefficients  and  of 
degree  n>l,  and  if  a  and  b  are  real  numbers  such  that  P(a) 
and  P(b)  have  opposite  signs,  then  the  polynomial  P(x)  has 
at  least  one  root  between  a  and  b.(22) 
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Example  A- 2 

Find  to  two  decimal  places  that  root  of  the  equation: 

s3  +  s2  -  12s  +  11  =  0  (A-7) 

which  lies  between  2  and  3. 

The  root  in  question  will  be  of  the  form  2.abc..., 
where  a,  b,  c,  ...  stand  for  its  decimal  figures. 
(1) :   Using  synthetic  division,  first  transform  equation  A-7 
into  an  equation  whose  roots  are  those  of  equation  A-7  dimin- 
ished by  2. 


1 

1 

-12 
6 

11 
-12 

11 

1 

3 
2 

-    6 
10 

-    1 

1 

5 

4 

2 

1   7 
The  transformed  equation  is: 

s3  +  7s2  +  4s  -  1  =  0  (A-8) 

Equation  A-8  will  have  the  root  . abc...  lying  between  0  and 
1. 

(2):   Substituting  the  values  s=  . 1,  .2,  .3,  ...  into  equa- 
tion A-8,  one  finds  that: 

f(.l)  <  0 
and 

f(.2)  >  0 
hence  equation  A-8  has  a  root  between  .1  and  .2,  so  that: 

a  =  1  (A-9) 
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(3):   Diminishing  the  roots  of  equation  A-8  by  .1  using  syn- 
thetic division: 

17     4      -1       .1 

.1    .71     .471 
I   771   3771   -  .529 

.1    .72 
1   772   57T3 

._1 

1   7.3 

The  transformed  equation  is: 

s3  +  7.3s2  +  5.43s  -  .529  =  0  (A-10) 

which  has  a  root  of  the  form  . Obc  .  .  . 

(4) :   Substituting  the  values  s  =  .01,  .02,  ...  into  equa- 
tion A-10,  one  finds  that: 

g(.08)  <  0 
and 

g(.09)  >  0   . 
Hence  equation  A-10  has  a  root  between  .08  and  ,09  so  that: 

b  =  8  . 
From  equation  A-9  and  A-ll,  the  root  of  equation  A-7 
which  lies  between  2  and  3  to  two  decimal  places  is: 

s  =  2.18  . 
The  process  could  be  continued  indefinitely,  thus  per- 
mitting one  to  compute  the  root  of  equation  A-7  between  2  and 
3  to  any  specified  number  of  decimal  places. 


101 


INITIAL  DISTRIBUTION  LIST 

No.  of  Copies 


1.  Defense  Documentation  Center  20 
Cameron  Station 

Alexandria,  Virginia   22314 

2.  Library  2 
Naval  Postgraduate  School 

Monterey,  California   93940 

3.  Director,  Naval  Communications  (OP-94)  1 
Department  of  the  Navy 

Washington,  D.  C.   20370 

4.  Professor  S.  R.  Parker  5 
Department  of  Electrical  Engineering 

Naval  Postgraduate  School 
Monterey,  California   93940 

5.  Lieutenant  T.  E.  Stone  1 
2  32  Anchorage  Road 

Annapolis,  Maryland   21401 


102 


UNCLASSIFIED 


Security  Classification 


DOCUMENT  CONTROL  DATA  -R&D 

itv  classification  of  title,   body  of  abstract  unci  indexing  annotation  must  be  entered  when  the  over. ill  report  is  cl.<?     ffied) 


1      ORIGINATING    ACTIVITY   (Corporate  uiithor) 


Naval  Postgraduate  School 
Monterey,  California   93940 


23.    REPORT    SECURITY     CLASSIMCA1ION 


Unclassified 


2b.    GROUP 


*      REPORT     TITLE 


Application  of  Associated  Polynomials  to  Circuits  and  Systems 


4     DESCRIPTIVE   NOTES  f7>pe  of  rfporl  <ind,inclusivf  dales) 

Thesis 


5     Au  THOR(S)  (First  name,  middle  initial,    last  name) 


STONE,  Thomas  Edward,  Lieutenant,  USN 


6       REPOR  T    DATE 


September    19  6  8 


la.     TOTAL    NO.    OF    PAGES 

103 


76.    NO.    OF    REFS 

23 


8a.    CONTRACT    OR    GRANT    NO. 


b.    PROJEC  T    NO. 


3a.    ORIGINATOR'S    REPORT    NUMBER(S) 


9b.   OTHER   REPORT   NOI5I  (Any  other  numbers   that  may  be  assigned 
this  report) 


10      DISTRIBUTION    STATEMENT 


mmmmim*BBbBta—SBa 


7iT£i£c\  rr.ii  gn 


■B 


OMHSBasasi 


93940 


II.    SUPPLEMENTARY    NOTES 


12.    SPONSORING    MILITARY    ACTIVITY 


Naval  Postgraduate  School 
Monterey,  California   93940 


1  3.  ABSTR  AC  T 


In  this  thesis,  the  theory  of  associated  oolynomials  is 
applied  to  linear  circuits  and  systems.   Associated  polyno- 
mials are  defined  and  used  in  the  presentation  of  a  computer 
algorithm  suitable  for  accomplishing  a  partial  fraction  ex- 
pansion including  repeated  roots.   Variations  of  the  coeffi- 
cients of  a  polynomial  are  related  to  the  variation  of  its 
roots  (and  vice-versa)  by  associated  polynomials.   These  re- 
sults are  used  in  a  root  solving  process  and  to  express  sen- 
sitivity coefficients  in  a  new  analytical  form.   Lastly,  a 
theory  of  compensating  parameter  adjustments  is  developed 
with  applications  to  self-adaptive  systems.   It  is  shown  that 
by  compensating  adjustments  of  any  two  parameters  in  the  sys- 
tem which  are  linearly  related  to  the  coefficients  of  the 
characteristic  equation,  adaptive  compensation  is  possible 
for  a  real  root  of  the  closed  loop  transfer  function  only 
and  that  for  complex  roots ,  the  real  and  imaginary  parts  may 
be  kept  invarient  but  not  both. 


DD  ,F°„R:.,1473 

S/N    01 01 -807-681  1 


(PAGE    1) 


UNCLASSIFIED 


103 


Security  Classification 


A-:ii408 


UNCLASSIFIED 


Security  Classification 


KEY     WO  R  DS 


LINK     A 


ROLE  W  T 


ASSOCIATED  POLYNOMIALS 
PARTIAL  FRACTION  EXPANSION 
SENSITIVITY  ANALYSIS 
SENSITIVITY  COEFFICIENTS 
SELF-ADAPTIVE  CONTROL  SYSTEMS 
VARIABLE  PARAMETERS 
CONTROL  SYSTEMS 


DD  ,FN°oR:,91473   back) 


S/N  0101-807-68?! 


104 


UNCLASSIFIED 


Security  Classification 


Ljaylorcl ; 

PAMPHLET  BINDER 

=S  Syracuse,  N.  Y. 
;  Stockton.  Calif,  j 


thesS764 


3  2768  0041 4845  2 


DUDLEY  KNOX  LIBRARY 


